arXiv:1501.02270v2 [astro-ph.GA] 21 Apr 2015 


Draft version April 22, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


CHAOS 1. DIRECT CHEMICAL ABUNDANCES FOR H II REGIONS IN NGC 628 

Danielle A. Berg^’^, Evan D. Skillman\ Kevin V. Croxall^, Richard W. Pogge^’'*, John Moustakas®, Mara 

Johnson-Groh® 

Draft version April 22, 2015 

ABSTRACT 

The CHemical Abundances of Spirals (CHAOS) project leverages the combined power of the Large 
Binocular Telescope (LBT) with the broad spectral range and sensitivity of the Multi Object Double 
Spectrograph (MODS) to measure “direct” abundances (based on observations of the temperature- 
sensitive auroral lines) in large samples of H II regions in spiral galaxies. We present LBT MODS 
observations of 62 H II regions in the nearby spiral galaxy NGC 628, with an unprecedentedly large 
number of auroral lines measurements (18 [O III] A4363, 29 [N II] A5755, 40 [S III]A63I2, and 40 [O II] 
AA7320, 7330 detections) in 45 H II regions. Comparing derived temperatures from multiple auroral 
line measurements, we find: (I) a strong correlation between temperatures based on [S III] A6312 
and [N II] A5755; and (2) large discrepancies for some temperatures based on [O II] AA7320, 7330 
and [O III] A4363. Both of these trends are consistent with other observations in the literature, 
yet, given the widespread use and acceptance of [O III] A4363 as a temperature determinant, the 
magnitude of the T[0 III] discrepancies still came as a surprise. Based on these results, we conduct a 
uniform abundance analysis prioritizing the temperatures derived from [S III] A63I2 and [N II] A5755, 
and report the gas-phase abundance gradients for NGC 628. Relative abundances of S/0, Ne/0, 
and Ar/0 are constant across the galaxy, consistent with no systematic change in the upper IMF 
over the sampled range in metallicity. These alpha-element ratios, along with N/0, all show small 
dispersions (cr ~ O.I dex) over 70% of the azimuthally averaged radius. We interpret these results as 
an indication that, at a given radius, the interstellar medium in NGC 628 is chemically well-mixed. 
Unlike the gradients in the nearly temperature-independent relative abundances, 0/H abundances 
have a larger intrinsic dispersion of ^ 0.165 dex. We posit that this dispersion represents an upper 
limit to the true dispersion in 0/H at a given radius and that some of that dispersion is due to 
systematic uncertainties arising from temperature measurements. 

Subject headings: galaxies: abundances - galaxies: spiral - galaxies: evolution 


1. INTRODUCTION 

The chemical evolution of galaxies, whereby succes¬ 
sive generations of stars enrich the interstellar medium 
(ISM) with products from stellar nucleosynthesis, is key 
to gaining a deeper understanding of galaxy evolution. 
H II regions can be used to study the absolute and rela¬ 
tive abundances in the ISM of spiral galaxies. Thus, spi¬ 
ral galaxies in the nearby universe, with low inclinations, 
provide the opportunity for studying and understanding 
the chemical evolution of galaxies. 

Chemical abundances have been widely studied in spi¬ 
ral galaxies using “strong-line” calibrations, where ra¬ 
tios of the strong forbidden lines are used as abundance 
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indicators. However, these measurements are only sta¬ 
tistical indications of chemical abundance, and are lim¬ 
ited by the parameter space of the calibrat i on sample 


(e.g. 

Ivan Zee & HavnesI 120061: IBresolinl 120071: lYin et al. 

2007 

;IBresolin et al.ll2009al:IStasihskall2010l:IAmorm et al. 

201C 

; iBerff et al.l ]2011[1. Because spiral galaxies ex- 


hibit abundance gradients and large ranges in temper¬ 
ature, density, and excitation space (among other pa¬ 
rameters), “strong-line” calibr ations are vulnerable to 
large systematic unce rtainties (jKewlev fc EllisonI l2008t 
iMoustakas et al.ll2010ll . In contrast, the “direct” method 
uses the temperature-sensitive ratio of auroral to forbid¬ 
den emission lines in order to determine the electron tem¬ 
perature and subsequent chemical abundances. 

Due to the inherent challenges of the direct method, 
detailed direct abundance studies exist for only a handful 
of spiral galaxie s (e.g., 20 H II regions w ith direct abun¬ 
dances in MlOl. lKennicutt et al.l (l^03bll: 61 H II regions 
with direct abundance in M33. IRosolowskv fc Sim^ 
(|2008f) ~). This relative dearth of large samples of di¬ 
rect abundances in the ISM of spiral galaxies prevents 
us from understanding key processes in the evolution of 
galaxies and limits our ability to fully understand: the 
chemical abundance gradients in spirals; the dispersion in 
the ISM abundances as a function of radius and environ¬ 
ment; the yields of elements from stellar nucleosynthesis 
with potential constraints on the stellar initial mass func¬ 
tion; the chemical evolution of galaxies at high redshift; 
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and thus, the chemical enrichment history of the uni¬ 
verse. In order to place these measurements on the same 
scale for comparison amongst galaxies, reliable and con- 
sistent derivations of abundances are needed (see, e.g., 
iMonstakas et al.ll201(m . 

With recognition that simply measuring a direct abun¬ 
dance is not the panacea to all H II region abundance 
measurement problems, assembling large, homogeneous 
datasets of electron temperature measurements repre¬ 
sent our best observational approach to putting all H II 
region abundance measurements on firmer ground. In¬ 
creasing the number of H II regions observed per galaxy 
allows for statistical approaches to questions about 
discrepant electron temperature measure ments (e.g., 
iKennicutt et al.l[2003bt iBinette et al.ll20I^. dispersions 
in abundances at a given ra.dius le.g.. IKennicutt et al.l 
1200, Sbl : iRosolowskv fc Simonll2008ll. and dispersions a s a 

function of environment le.g.. IKennicutt et al.ll200.‘Ib[l . 

1.1. The CHAOS Project 

The CHemical Abundances Of Spirals (CHAOS) 
projectQ was undertaken in order to build a large 
database of high quality H II region spectra from nearby 
spiral galaxies using the Multi-Ob ject Double Spectro¬ 
graphs IMODS: IPogge et al.ll20I0ll on the Large Binoc¬ 
ular Telescope (LBT). The MODS have been designed 
to optimally obtain high quality spectra across fields of 
view (6'x6') comparable to the extents of nearby spiral 
galaxies. Furthermore, the LBT and MODS combina¬ 
tion provides the balance between sensitivity, resolution, 
and wavelength coverage necessary to measure all emis¬ 
sion lines relevant to oxygen, nitrogen, and alpha-element 
abundance determinations and more. In particular, the 
optical design of MODS has been optimized for high 
throughput in the red and blue channels to produce high- 
quality spectrophotometry from the atmospheric UV cut¬ 
off (^340 nm) to the silicon detector cutoff at ^1/rm. 

The Syitzer Infrared Nearby Galaxies Survey (SINGS; 
IKennicutt et al.l l2003all galaxies are arguably the best 
understood spiral galaxies outside of the Local Group 
due to having the most complete multi-wavelength ob¬ 
servations available. Due to the Spitzer observations and 
large ancillary data sets, these galaxies have resolved 3.6- 
160 pLiR imaging, 5-40 ^m IRS spectroscopy of the central 
regions plus select extra-nuclear H II regions, and high- 
quality H I and GO gas maps with detailed 2-D rotation 
curves. Additionally, observations by Herschel under the 
Key Insights on Nearby Galaxies: a F ar- Infrared Sur¬ 
vey w ith Herschel project (KINGFISH; IKennicutt et al.l 
I20III1 provide PACS and SPIRE imaging, as well as far- 
IR spectroscopy of selected H II regions. We defined 
the GHAOS sample by identifying the optimal candi¬ 
dates within the SINGS sample: we restrict our sam¬ 
ple to nearly face-on, luminous (Mb < —18 mag), spiral 
galaxies (Hubble type T > 0) accessible from the north¬ 
ern hemisphere with the LBT (declination > —5°), with 
minimal transit airmasses (< 1.1). For our final target 
sample, 13 high priority SINGS galaxies were selected. 

The collection of a large sample of H II regions is driven 
by the need to precisely measure the form of the ra¬ 
dial gradient, including deviations from simple linear fits, 
and the dispersion in abundance at fixed radius. Ideally, 
roughly 10 H II regions in approximately 10 independent 
radial bins, or ^ 100 H II regions per galaxy, would be 


required to achieve the desired precision. Given the steep 
H II region luminosity function, and the non-random dis¬ 
tribution of H II regions in a spiral galaxy, this ideal sta¬ 
tistical sample is not obtainable for all galaxies, but rep¬ 
resents a reasonable goal. At present, multiple fields have 
been observed in 9 of the 13 galaxies, already more than 
doubling the number of H II regions in spiral galaxies 
with sufficient quality spectra to provide accurate mea¬ 
surements of the physical conditions (i.e., temperatures, 
densities, and ionization parameters) and absolute and 
relative chemical abundances. 


1.2. NGC 628 Properties 

We present LBT observations of H II regions in the 
first CHAOS target, NGC 628. NGC 628 (M 74) 
is a late-type giant Sc spiral galaxy with a system¬ 
atic velocity of 656 km s~^. W e adopt a distance 
of 7.2 Mdc (I Van Dvk et al.l I2006IFI and an inclination 
of z « 5° (iShostak fc van der KruitI Il984ll . with a re¬ 
sulting scale of 35 pc arcsec”^. The adopted proper¬ 
ties for NGC 628 are provided in Table [fl NGC 628 
is an excellent target due to its small inclination, ex¬ 
tended structure, and undisturbed optical profile. The 
gas-phase oxygen abundance of NGC 628 has been 
studied previ o usly using long-slit or fiber spectroscopy 
(iTalentI 119831: ^Call et all 19851: iZaritskv et al.l 11 994 


Ferguson et al. 19981: ^n Zee et al.ll998bl: Bresolin et al l 


19991: Castellanos et al. 20021: iMoustal^s et al. 20101 : 


Gusev et al.l 120121 : iCedrds et al.l 1^121: Berg et al. 1201311. 

and integral field spec troscopy (|Rosales-Ortega et al l 
I2011l : [&oxall et aLll2013ll . All together, a combined total 
of 18 distinct H II region direct abunda nce measurements 
exist in the literature for NGC 628 (ICastellanos et aO 
120021 : IBerg et HI 120131: ICroxall et al.l l2013rFI . The 
CHAOS observations of NGC 628 present a three-fold 
increase in the total number of direct H II region abun¬ 
dances and allow, for the first time, the spatial variations 
and dispersions in absolute and relative ISM abundances 
to be determined from the direct method. 

This paper is organized as follows: In Section [2] we de¬ 
scribe the LBT observations for NGC 628, the spectral 
data processing, the emission line measurements, and the 
reddening corrections. In Section [3] we discuss the deriva¬ 
tion of the temperatures and densities of the H II regions, 
including the choice of atomic data and comparison of 
the derived temperatures. Our derived abundances are 
presented in Section Ul where we introduce our uniform 
abundance analysis method and report ionic and abso¬ 
lute abundances for NGC 628. We further interpret the 
CHAOS abundance gradients for NGC 62 8 by consider¬ 
ing electron temperature discrepancies (§ EM , co mpar - 
ing to the abundance gradients in the literature (§ EM, 
and assessing the physical basis for dispersions at a given 
galactocentric radius. Finally, we summarize our find¬ 
ings in Section ini We also present additional analyses in 
two Appendices: The direct oxygen abundance gradient 
based on [O III] A4363 electron temperature measure¬ 
ments is presented in Appendix Indicators of temper- 
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TABLE 1 

Properties of NGC 628 


Property 

NGC 628 

R.A. 

01:36:41.747 

Dec. 

15:47:01.18 

Type 

Scl 

Adopted D (Mpc) 

7.2±1.0l 

ms (mag) 

9.952 

Redshift 

0.002192 

Inclination (degrees) 

53 

P.A. (degrees) 

124 

i ?25 (arcmin) 

5.25® 

R 25 (kpc) 

10.95 


Note. — Adopted properties 
for NGC 628. Rows 1 and 2 
give the RA and Dec of the opti¬ 
cal center in units of hours, min¬ 
utes, and seconds, and decrees, ar- 
cminutes, and arcseconds respec¬ 
tively. Row 5 lists redshifts taken 
from the NASA/IPAC Extragalactic 
Database. Row 8 gives the optical 
radius at the B 25 mag arcsec”^ of 
the system. Row 9 gives the op¬ 
tical radius of the galaxy given the 
adopted distance. R efere nces: (1) 

I Van Dvk et al.l 11200611: (2) ILee et al. 

^0^1); (3) [Shosta^^ van der^ruin 

('4') lEeus^^^hT 1120091): (51 
iKeJail Jhnm 

ature discrepancies are presented in Appendix [BJ 

2. CHAOS SPECTROSCOPIC OBSERVATIONS 
2.1. Observation Plan 

Optical spectra of NGC 628 were acquired with 
MODSl on the LBT as part of the CHAOS project on 
the UT dates of 2012 October 14 and 16. The pri¬ 
mary goal was to obtain high signal-to-noise spectra, 
with detections of intrinsically faint auroral lines (e.g., 
[O III] A4363, [N II] A5755, JS II] A6312) at a signifi¬ 
cance of 3(7 or higher. The multi-object mode of MODS, 
which uses custom-designed, laser-milled slit masks, al¬ 
lows spectra of many H II regions to be obtained simul¬ 
taneously. Broad-band and Ha continuum-subtracted 
SINGS images for NGC 628 (|Muhoz-Mateos et al.ll2009ll 
were used to identify target H II regions, as well as align¬ 
ment stars, and determine accurate astrometry for the 
masks. H II regions were selected to achieve large radial 
coverage of the optical disk, prioritizing knots of high¬ 
est Ha surface brightness. We observed three multi-slit 
masks, each containing ^20 slits, which spanned the ra¬ 
dial and azimuthal extent of the optical disk of NGC 628. 
All H II region slits are 1.0" wide, but lengths vary be¬ 
tween 8 — 20" depending on the size of the targeted H II 
region and proximity of other slits on the mask. Blue 
and red spectra were obtained simultaneously using the 
G400L (400 lines mm"!, R~1850) and G670L (250 lines 
mm“^, R~2300) gratings, respectively. This combina¬ 
tion results in a broad spectral coverage extending from 

Funding awarded by NSF Grant No. 1109066. 

^ IVan Dvk et aO II2006I1 reports a value from supernovae distance 
determinations. Current value from the literature for the distance 
of NGC 628 range from 7.2 — 9.95 Mpc. 

® 25 observations may be totaled from these three studies, how¬ 
ever, several originate from the same H II region complexes. Thus, 
only 18 non-duplicate region observations exist. 



Fig. 1.^— Continuum-subt racted Ho SINGS image of NGC 628 
IIMuhoz-Mateos et all [20091) . The footprints of CHAOS slits are 
overlaid, where red, green, and blue respectively represent the Field 
1, 2, and 3 slit positions observed at the LBT. The slit positions 
targeted H II regions, although not always centered in order to 
maximize effective usage of mask real estate. See Table [2] for more 
details. 

3200 - 10,000 A, that was linearly re-sampled to 0.5 A 
per pixel and full width half maximum resolution of ^ 2 

A. 

Each mask field was observed for six 1200 second ex¬ 
posures, for a total integration time of two hours. The 
masks were designed with slits at a fixed position angle 
which approximated the parallactic angle at the mid¬ 
point of the longest possible observation window in a 
night. For NGC 628, during the observing run of UTC 
2012 October 2—17, we selected the parallactic angle of 
PA = —235.2° for Fields 1 and 3, and PA = —55.1° for 
Field 2, which is just a rotation of 180° to allow for the 
use of different guide stars and thus a shift of the field 
center. This, in addition to observing the galaxies at air- 
masses less than 1.5, served to minimize the wavelength- 
depende nt light loss due to differential atmospheric re¬ 
fraction (|FiliDDenkcilll982f) . On the nights that NGC 628 
was observed, the sky conditions were optimal: clear, 
low wind, and low humidity. Additionally, we aimed for 
approximately arcsecond seeing. Fields 1, 2, and 3 had 
average seeing values over their 2 hour observation win¬ 
dows of 1.1, 1.25, and 0.85 arcseconds, respectively. 

Figure [T] shows t he Ha continuum-subtract ed SINGS 
image for NGC 628 (|Munoz-Mateos et al.ll200^ . with slit 
positions from Table [5] shown as colored lines for each 
mask field. Throughout this work, all target locations 
(extraction of H II region; not slit center) are referred to 
by their offsets in right ascension and declination from 
the center of NGC 628, which is taken to be 01^36"‘41.7'* 
-|-15°47'01". Slit locations and offsets for all observed 
H II regions are listed in Table [H 

2.2. Spectral Reductions 
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TABLE 2 

Logs for NGC 628 LBT Observations 


H II 

Region 

R.A. 

(2000) 

Dec. 

(2000) 

Rg. 

(arcsec) 

R/R.25 

Rg 

(kpc) 

Offset (R.A., Dec.) 
(arcsec) 

Auroral Line Detections 

[0 iiij-^Tn]—[STo]—[DTIT 

(y iij 

Total Detections: 

18 

29 

40 

40 

0 

NGC628-12.7+2.1 

1:36:40.92 

15:47:02.06 

12.037 

0.038 

0.42 

-12.7,+2.1 






NGC628+38.6-17.8 

1:36:44.47 

15:46:42.19 

43.863 

0.139 

1.53 

+38.6,-17.8 






NGC628-48.8+0.7 

1:36:38.42 

15:47:00.72 

48.194 

0.153 

1.68 

-48.8,+0.7 






NGC628+38.9-26.7 

1:36:44.49 

15:46:33.32 

48.628 

0.154 

1.69 

+38.9,-26.7 






NGC628-44.7+27.3 

1:36:38.70 

15:47:27.34 

51.311 

0.163 

1.79 

-44.7,+27.3 






NGC628-33.2+45.8 

1:36:39.50 

15:47:45.80 

55.297 

0.176 

1.93 

-33.2,+45.8 






NGC628-30.7-46.5 

1:36:39.67 

15:46:13.48 

56.337 

0.179 

1.97 

-30.7,-46.5 






NGC628+48.0-32.0 

1:36:45.13 

15:46:27.98 

59.216 

0.188 

2.07 

+48.0,-32.0 






NGC628-35.9+57.7 

1:36:39.31 

15:47:57.72 

66.695 

0.212 

2.33 

-35.9,+57.7 


/ 

/ 



NGC628-55.5-40.4 

1:36:37.96 

15:46:19.64 

68.798 

0.218 

2.40 

-55.5,-40.4 






NGC628+49.8+48.7 

1:36:45.25 

15:47:48.72 

69.474 

0.221 

2.43 

+49.8,+48.7 


/ 

/ 

/ 


NGC628-8.3-72.8 

1:36:41.23 

15:45:47.18 

74.386 

0.236 

2.59 

-8.3,-72.8 






NGC628-73.1-27.3 

1:36:36.73 

15:46:32.68 

77.973 

0.248 

2.72 

-73.1,-27.3 



/ 

/ 


NGC628-76.2+22.9 

1:36:36.52 

15:47:22.91 

78.816 

0.250 

2.75 

-76.2,+22.9 


/ 

/ 

/ 


NGC628+18.5+79.0 

1:36:43.08 

15:48:19.02 

80.178 

0.255 

2.80 

+18.5,+79.0 






NGC628-36.8-73.4 

1:36:39.25 

15:45:46.65 

82.827 

0.263 

2.89 

-36.8,-73.4 


/ 


/ 


NGC628-73.1-44.8 

1:36:36.74 

15:46:15.18 

85.853 

0.273 

3.00 

-73.1,-44.8 






NGC628+68.5+53.4 

1:36:46.54 

15:47:53.35 

86.834 

0.276 

3.04 

+68.5,+53.4 


/ 

/ 

/ 


NGC628-87.5-10.8 

1:36:35.74 

15:46:49.16 

87.870 

0.279 

3.07 

-87.5,-10.8 






NGC628+81.6-32.3 

1:36:47.45 

15:46:27.69 

89.257 

0.283 

3.11 

+81.6,-32.3 



/ 



NGC628-68.5+61.7 

1:36:37.05 

15:48:01.74 

91.116 

0.289 

3.18 

-68.5,+61.7 


/ 

/ 

/ 


NGC628+76.9-49.6 

1:36:47.13 

15:46:10.42 

93.074 

0.295 

3.25 

+76.9,-49.6 


/ 

/ 

/ 


NGC628+94.2+4.0 

1:36:48.33 

15:47:03.99 

95.361 

0.303 

3.33 

+94.2,+4.0 






NGC628+75.4+66.0 

1:36:47.02 

15:48:06.01 

100.177 

0.318 

3.50 

+75.4,+66.0 






NGC628-13.1+107.5 

1:36:40.89 

15:48:47.50 

107.076 

0.340 

3.73 

-13.1,+107.5 


/ 

/ 

/ 


NGC628+53.5-104.0 

1:36:45.51 

15:45:16.04 

118.496 

0.376 

4.13 

+53.5,-104.0 


/ 


/ 


NGC628-35.7+119.6 

1:36:39.33 

15:48:59.59 

123.566 

0.392 

4.31 

-35.7,+119.6 


/ 

/ 

/ 


NGC628-20.3+124.6 

1:36:40.40 

15:49:04.62 

125.033 

0.397 

4.36 

-20.3,+124.6 


/ 

/ 

/ 


NGC628-59.6-111.6 

1:36:37.67 

15:45:08.43 

127.196 

0.404 

4.44 

-59.6,-111.6 


/ 

/ 

/ 


NGC628+61.2+113.5 

1:36:46.04 

15:48:53.45 

128.290 

0.407 

4.48 

+61.2,+113.5 


/ 

/ 

/ 


NGC628-128.4+13.1 

1:36:32.90 

15:47:13.09 

128.700 

0.409 

4.50 

-128.4,+13.1 






NGC628+42.6-120.7 

1:36:44.75 

15:44:59.30 

129.490 

0.411 

4.52 

+42.6,-120.7 


/ 

/ 

/ 


NGC628+131.9+18.5 

1:36:50.94 

15:47:18.53 

134.275 

0.426 

4.69 

+131.9,+18.5 

/ 

/ 

/ 

/ 


NGC628+71.2+121.3 

1:36:46.74 

15:49:01.32 

140.117 

0.445 

4.90 

+71.2,+121.3 






NGC628+125.4-62.4 

1:36:50.49 

15:45:57.61 

141.767 

0.450 

4.95 

+125.4,-62.4 



/ 

/ 


NGC628-130.9+71.8 

1:36:32.73 

15:48:11.80 

148.566 

0.472 

5.18 

-130.9,+71.8 


/ 

/ 

/ 


NGC628+131.7-70.2 

1:36:50.93 

15:45:49.80 

151.024 

0.479 

5.27 

+131.7,-70.2 



/ 

/ 


NGC628+151.0+22.3 

1:36:52.26 

15:47:22.33 

153.790 

0.488 

5.37 

+151.0,+22.3 

/ 

/ 

/ 

/ 


NGC628-157.9-0.3 

1:36:30.86 

15:46:59.71 

157.687 

0.501 

5.51 

-157.9,-0.3 


/ 

/ 

/ 


NGC628-24.5-155.6 

1:36:40.10 

15:44:24.39 

158.586 

0.503 

5.54 

-24.5,-155.6 


/ 

/ 

/ 


NGC628-129.8+94.7 

1:36:32.81 

15:48:34.73 

159.867 

0.508 

5.58 

-129.8,+94.7 


/ 

/ 

/ 


NGC628+158.5+9.2 

1:36:52.78 

15:47:09.20 

160.025 

0.508 

5.58 

+158.5,+9.2 




/ 


NGC628+78.6-137.6 

1:36:47.24 

15:44:42.37 

160.143 

0.508 

5.59 

+78.6,-137.6 




/ 


NGC628+140.3+82.0 

1:36:51.52 

15:48:21.99 

162.919 

0.517 

5.68 

+140.3,+82.0 



/ 



NGC628-42.8-158.2 

1:36:38.84 

15:44:21.79 

164.834 

0.523 

5.75 

-42.8,-158.2 

/ 

/ 

/ 

/ 


NGC628+147.9-71.8 

1:36:52.05 

15:45:48.18 

166.229 

0.528 

5.80 

+147.9,-71.8 


/ 


/ 


NGC628+163.5+64.4 

1:36:53.13 

15:48:04.34 

176.489 

0.560 

6.16 

+163.5,+64.4 


/ 

/ 

/ 


NGC628-4.5+185.6 

1:36:41.49 

15:50:05.63 

184.519 

0.586 

6.44 

-4.5,+185.6 


/ 

/ 

/ 


NGC628+176.7-50.0 

1:36:54.04 

15:46:09.96 

185.425 

0.589 

6.47 

+176.7,-50.0 



/ 

/ 


NGC628-76.2-171.8 

1:36:36.52 

15:44:08.22 

188.715 

0.599 

6.59 

-76.2,-171.8 

/ 


/ 

/ 


NGC628+31.6-191.1 

1:36:43.99 

15:43:48.85 

195.131 

0.619 

6.81 

+31.6,-191.1 

/ 

/ 

/ 

/ 


NGC628-200.6-4.2 

1:36:27.90 

15:46:55.81 

200.634 

0.637 

7.00 

-200.6,-4.2 

/ 

/ 

/ 

/ 


NGC628-184.7+83.4 

1:36:29.00 

15:48:23.36 

202.182 

0.642 

7.06 

-184.7,+83.4 

/ 

/ 

/ 

/ 


NGC628-206.5-25.7 

1:36:27.49 

15:46:34.31 

208.222 

0.661 

7.27 

-206.5,-25.7 

/ 


/ 

/ 


NGC628-90.1+190.2 

1:36:35.56 

15:50:10.17 

209.327 

0.665 

7.31 

-90.1,+190.2 

/ 

/ 

/ 

/ 


NGC628-168.2+150.8 

1:36:30.14 

15:49:30.77 

225.204 

0.715 

7.86 

-168.2,+150.8 

/ 


/ 

/ 


NGC628+232.7+6.6 

1:36:57.92 

15:47:06.58 

234.411 

0.744 

8.19 

+232.7,+6.6 

/ 


/ 

/ 


NGC628+237.6+3.0 

1:36:58.26 

15:47:02.93 

239.203 

0.759 

8.35 

+237.6,+3.0 

/ 


/ 

/ 


NGC628+254.3-42.8 

1:36:59.42 

15:46:17.11 

259.853 

0.825 

9.07 

+254.3,-42.8 

/ 


/ 



NGC628+252.1-92.1 

1:36:59.26 

15:45:27.85 

270.527 

0.859 

9.45 

+252.1,-92.1 

/ 

/ 

/ 



NGC628+261.9-99.7 

1:36:59.94 

15:45:20.24 

282.451 

0.897 

9.89 

+261.9,-99.7 

/ 





NGC628+265.2-102.2 

1:37:00.17 

15:45:17.76 

286.416 

0.909 

10.00 

+265.2,-102.2 

/ 





NGC628+289.9-17.4 

1:37:01.89 

15:46:42.55 

292.382 

0.928 

10.20 

+289.9,-17.4 

/ 


/ 

/ 


NGC628+298.4+12.3 

1:37:02.47 

15:47:12.26 

300.442 

0.954 

10.49 

+298.4,+12.3 

/ 


/ 

/ 



Note. — Observing logs for H li regions observed in NGC 628 using MODS of the LBT on the UT dates of 14 and 16 October 2012. Each field 
was observed over an integrated exposure time of 1200 s on clear nights, with on average ~ l" seeing, and airmasses less than 1.5. H ll region IDs 
are listed in Column 1, labeled according to their offsets, and ordered by increasing galactocentric radius. The right ascension and declination of 
the individual H li regions are given in units of hours, minutes, and seconds, and decrees, arcminutes, and arcseconds respectively. The H ll region 
distances from the center of the galaxy in arcseconds, fraction of R 2 b, and in kpc are listed in the Columns 4-6. Additionally, the offsets from galaxy 
center R.A. and Dec. are given in arcseconds in Column 7. Finally, Columns 8-12 highlight which regions have [O III] A4363, [N il] A5755, [S ill] 
A6312, O ii] AA7320, 7330, and S ii] AA4068, 4076 auroral lines detections at the Aa significance level. 
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Spectra are reduced and analyzed using the beta- 
release version of the MODS reduction pipeline^. 
Bias subtraction and flat-fielding were performed using 
python scripts designed to accommodate the interlaced 
data read-out. The pixel flat is assembled from slit¬ 
less lamp flats of the incandescent continuum lamp. Slit 
definitions, wavelength calibration, and creation of vari¬ 
ance images are performed using new tools written to 
work within the XIDiO reduction package. Slit loca¬ 
tions on two-dimensional images are defined by mapping 
laser-milled mask positions and wavelength to pixel co¬ 
ordinates on the CCD. Wavelength calibrations are very 
stable, modulo a small flexure shift of order a few pix¬ 
els. Spectra of Hg -|- Ar, Ne, and Xe -I- Kr lamps are 
taken through the slit masks during a calibration se¬ 
quence performed at the beginning of the observing run. 
Because the MODS spectrographs are remarkably sta¬ 
ble, the optical path is mapped by using the images of 
lamps through specially designed pinhole masks. With 
this information, we generate a transformation from an 
X/Y location on the mask to a pixel on the detector as 
a function of wavelength, and this geometric transforma¬ 
tion allows the construction of an initial trace. The initial 
spectral traces are corrected for atmospheric refraction 
using the hour angle and airmass corresponding to the 
mid-point of the observation. This extraction technique 
was designed for and is effective for use with very faint 
spectra where continua are faint or not visible. No slit- 
losses are observed above 3500 A and so no correction is 
made to this end. 

Due to the large angular size of NGC 628, the diffuse 
ionized medium fills the majority of each field of view. 
In order to facilitate accurate sky subtraction, additional 
“sky” slits are placed well away from both H II regions 
and the main body of the galaxy. Due to the off-axis na¬ 
ture of the sky-slits, they are more affected by aberration 
than slits targeting H II regions. A sky frame was created 
by fitting a two-dimensional B-spline to the background 
sky in the two-dimensional spectra, i.e., each slit is fit in¬ 
dependently with a local sky. The continuum flux from 
the sky-slit is scaled to the background continuum mea¬ 
sured locally in each slit. Subsequently, the background 
under the strong nebular lines is interpolated from the 
scaled sky-slit to create the final sky-model that is sub¬ 
tracted. 

One-dimensional spectra are corrected for atmospheric 
extinction and flux calibrated based on obse rvations of 
CALSPEC flux standard stars (lBohlinir2010ll . The flux 
standards are observed on each night science data are 
obtained. We confirm the high level of fidelity across the 
spectra by comparing the continuum levels of the red 
and blue spectra and find that corrections for matching 
the two sides are on the order of ^ 1%. An example 
flux-calibrated spectrum is shown in Figure [2] 

2.3. Emission Line Measurements 

Emission line strengths are measured using various ex¬ 
isting tools described in iCroxall et al.l (1201311 . The un- 

The MODS reduction pipeline was developed by 
Kevin Croxall with funding from NSF Grant AST-1108693. 
Details can be found at - http://www.astronomy.ohio- 
state.edu/MODS/Software/modsIDL /. 

http://www.ucolick.org/xavier/IDL/. 


derlying continuum of our MODSl spectra are mod- 
eled using the STARL IGHlP^ s pectral synthes i s code 
(|Fernandes et al.l 1200511 . We fit iBruzual et al.l l)2013H 
models to the blue continuum red-ward of the [O II] 
A3727 line, after masking n ebular emission lines and 
Wolf-Rayet features, using a iCalzetti et al.l (j2000ll red¬ 
dening law. Subsequently, emission lines are fit simulta¬ 
neously with stellar and nebular continuum components. 
When automatically fitting emission lines, we initially as¬ 
sume these lines are well described by a Gaussian profile. 
Once continuum levels and line locations are measured, 
the final emission line values are computed by integrating 
the flux over each profile width. This step is important 
with respect to the high resolution of MODS data, be¬ 
cause nebular structure can be seen in the emission line 
profiles such that they are not well-characterized by a 
Gaussian function. To ensure that noise spikes and mul¬ 
tiple features within the same line are not mis-identified 
and fit, emission lines are required to have the same bulk 
kinematics and FWHM. We validate the automated line 
fits by measuring integrated fluxes by hand using the 
SPLOT routine within IR APF^ for all 62 extracted one¬ 
dimensional spectra. This exercise confirms the reliabil¬ 
ity of the automated method. 

In the past, standard errors associated with a given 
line strength often foll owed a relatively sim ple, general¬ 
ized prescription le.g.. iSkillman et aI1il994f) . Here, the 
error associated with each line detection is a combina¬ 
tion of the measured variance extracted from the two- 
dimensional variance image, the error associated with the 
sensitivity function, the Poisson noise in the continuum, 
the read noise of the detector, sky noise, a flat field¬ 
ing calibration error, the error in the continuum place¬ 
ment, and uncertainty associated with reddening. For 
weak lines, the uncertainty is dominated by error from 
the continuum subtraction, whereas the lines with flux 
measurements much stronger than the rms noise of the 
continuum have errors dominated by error in flux cali¬ 
bration and reddening. To account for the inherent un¬ 
certainty in t he flux cal ibration based on standard star 
observations (jOkdllQQflll . a minimum uncertainty of 2% 
is added in quadrature to the emission line flux uncer¬ 
tainties. Note that this is a conservative estimate as the 
use of HST has allowed improved ca librations of flux 
standards in GALSPEG (iBohlinI [201011 . The minimum 
2% emission line uncertainty forces a floor of ^ 50 K 
uncertainty for the electron temperature calculations. 

In this study, we focus on the most reliable spectra 
wherein at least one temperature-sensitive auroral line 
or doublet (i.e., [O III] A4363, [N II] A5755, [S III] A6312, 
or [O II] AA7320,7330) is detected at a line strength 
> 3(7. Note that this criterion is enforced to ensure that, 
despite best spectral reduction practices, contamination 
from other lines, continuum noise, remnant dichroic ef¬ 
fects, etc. do not lead to false auroral line detections. In 
addition to this precaution, each auroral line is examined 
by eye to ensure no valid detections are excluded and no 
suspicious detections pass the cut. By this criterion, out 
of 62 observed H II regions, we detect at least one auroral 
line in 45 spectra. 

http://astro.ufsc.br/starlight/. 

IRAF is distributed by the National Optical Astronomical 
Observatories. 
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NGC 628-42.8-158.2 





Wavelength (A) 


Fig. 2.— Sample LET spectrum of an H II region in NGC 628 with auroral line detections at a strength of 3cr or greater. The spectrum 
is broken up into three sections: The top panel is the blue region, where the high quality of the spectrum is evident by the high signal 
to noise Balmer series (highlighted in the expanded box). The middle panel corresponds to the middle portion of the spectrum and the 
bottom panel to the red region of the spectrum. The strength of the auroral lines used in this paper [O III] A4636, [N II] A5755, [S III] 
A6312, and [O II] AA7320, 7330) are highlighted in the corresponding expanded boxes. 
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2.4. Reddening Corrections 

The relative intensities of the Balmer lines are nearly 
independent of both density and temperature, so they 
can be used to solve for the reddening. The spec¬ 
tra are de-reddened using an iterat ive application of 
the reddening law of lO’Donnelll (|1994f) . parametrized by 
Av = 3.1 E{B — V). An initial reddening estimate is de¬ 
termined from the Balmer line ratios assuming a temper¬ 
ature of Te = 10^ K and a density of rie = 10^ cm“^. Us¬ 
ing this reddening, we determine an initial estimate of the 
electron temperature, which is an error weighted average 
of the electron temperatures calculated from the auroral 
lines. This new electron temperature is then used itera¬ 
tively to modify the input assumptions for the reddening 
code until the average electron temperature changes by 
less than 20 K. The final reddening estimate is an error 
weighted average of the individual reddening values de¬ 
termined from the H5/H/3, Hy/H/S, and Ha/H/? ratios. 
As the lines are fit simultaneously with stellar popula¬ 
tion models, further corrections for underlying Balmer 
absorption are unnecessary. 

Following iLee fc SkillmanI (j2n04[l . the reddening value 
can be converted to the logarithmic extinction at iJ/3 as 
C(H/3) = lA‘iE{B — V). Because of its high Galactic 
latitude, the foreground extinction for NGC 628 is ex- 
pecte d to be low, with an approxi mate value of G(H/3) = 
0.10 (jSchlaflv fc Finkbeinerl[2(TTTh . For NGG 628 we find 
a range in C(H/3) of ~ 0.0 to 0.6. Because the MODS 
spectra extend out to the 1 /rm cutoff, we observe the 
Paschen series and so perform an independent calculation 
of the reddening parameter which confirms the Balmer 
determinations. 

Reddening corrected line intensities measured from the 
observed H II regions in the target fields are reported in 
TableO (the full table is available online). Various optical 
emission line ratios can be used to distinguish emission 
regions which are ex cited by different mechanisms (e.g., 
[Baldwin et al.lFToSlf ) . All 45 regions in the CHAOS sam¬ 
ple clearly lie within the typical ranges observed for H II 
regions, ruling out contributions from active galactic nu¬ 
clei and supernovae remnants, and allowing the standard 
direct analysis of physical conditions in H II regions. 

3. PHYSICAL CONDITIONS 
3.1. Three-Zone Temperature Determinations 

A simple H II region can be modeled by three sepa¬ 
rate volumes: a low-, intermediate-, and high-ionization 
zone. Accurate H II region abundance determinations 
require reliable electron temperature measurements for 
each volume. This is typically done by observing a 
temperature-sensitive auroral-to-strong-line ratio. The 
[O III] I(AA4959,5007)/I(A4363) ratio is expected to re¬ 
flect the temperature in the high ionization zone. Sim¬ 
ilarly, the [N II] I(AA6548,6484)/I(A5755) and the [O II] 
I(AA7320,7330)/I(A3726,3728) ratios are expected to re¬ 
flect the temperature in the low ionization zone. S and 
Ar ions do not neatly separate into these two zones and 
thus require an intermediate ionization zone, which can 
be measured by the (S III] I(AA9069,9532)/I(A6312) ra¬ 
tio. Since the redder [S III] lines can be affected by at¬ 
mospheric absorption, we employ the theoretical ratio of 
A9532/A9069 = 2.51 to make an upward correction to 
the weaker of the two lines if a deviation of more than 


2% from the theoretical ratio is observed. 

Following the precedent of iBresolin et al.l (|2009af) . we 
have adopted the current best set of atomic data from 
the literature as reported in Table IH For the H II re¬ 
gions with detected auroral lines, we determine electron 
temperatures using the reddening corrected line-ratios 
and appropriate atomic data, assuming the ions are well- 
approximated by a 5-level atoirP^. Once a direct elec¬ 
tron temperature is determined, the physical conditions 
of the other zones are needed to com plete the H II re¬ 
gion picture. iGarnettI (hereafter G92: 1199^ relates the 
direct temperatures of different ionization zones using 
photoionization models: 

T[S III] = 0.83 X T(0 III] -L 1700 K (1) 

T(N ii] = 0.70 X T(0 III] -L 3000 K. (2) 

These relationships are valid for temperatures typical of 
H II regions: Tg = 2,000—18,OOOK. Observational ev¬ 
idence supporting these rel ationships can be foun d in 
the literature. For exam ple. IBresolin et al.l (j2009a[) and 
iKennicutt et al.l (|2003bf l find a tight correlation between 
H II region observations and the G92 T[S III ] versus 
T[0 III] relationship and lEsteban et al.l (|2009f) find a 
high correlation coefficient for the T[N II] versus T]0 III] 
G92 model. With the assumption that T[N II] = T[0 II], 
since N"*" and O'*" have comparable ionization energies, 
we now have the tools to determine temperatures in all 
three ionization zone for all ions. 


3.2. Impact Of Updated Atomic Data 

Recent updates to atomic data are of particular 
importance to this paper as changes in the Ein¬ 
stein A-coefficients for [S III] result in a system¬ 
atic sh ift of T[S III] to coo ler temperatures. In 
MIDI, IKennicutt et al.l (j2003bL hereafter K03) mea¬ 
sured 16 regions with both T[0 III] and T[S III] 
and found a tight correlation around relationship de¬ 
fined by_G92^^£^ used collision strengths for S++ 
from iTaval fc Guotal (Il999f) and Einstein A coefficients 
from Kauffman fc Sugar ~Tl986lf . Contemporaneously, 
iPerez-Montero fc Dfad (120031 ) noted a T[S III]-T[0 III] 
theoretical relationship that differs from the standard 
G92 relationship, where the d eparture is mos t ly du e 
to introducing atomic data from iTaval fc Guota (1199911. 
This is also a factor in the analysis bv iBinette et af 


(2012), wh o use the abundan c es fro m iHagele et al 


(2006) and IPerez-Montero et al.l (j2006l ). both of whom 
emplo y the [S III] collision strengths from iTaval fc Guiit^ 
(|1999D . Since then, further derivations of atomic coeffi- 
ci ents for have been pu blished, such as the work 
of iFroese Fischer et al.l (l 2006f). whose Einstein A-values 
we adopt. iMendoza fc Bautistal (12014( 1 performed exten¬ 
sive computations of new collisional data sets for [O III], 
[O II], [N II], ]S III], [S II], and [Ar III] and compared to 
other data sets of effective collision strengths computed 
within the past two decades. They found that statis¬ 
tical consistency amongst data sets is typically around 
^ 20 — 30%, where the inherent dispersion in forbid¬ 
den transitions is due to the strong sensitivity of the the 
resonance contribution to the scattering numerical ap¬ 
proach, target atomic model, energy mesh, and small 


https: / / github.com / moustakas/impro 
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TABLE 3 

Emission-Line Intensities and Equivalent Widths for LBT Observations of H II regions in NGC 628 


iw/nm 


Ion 

H14 A3721- 

[O ii] A3726 
[O ii] A3728 
H13 A3734 
H12 A3750 
Hll A3770 
HIO A3797 
He I A3819 
H9 A3835 
[Ne III] A3868 
He I A3888 
H8 A3889 
He I A3964 
[Ne III] A3967 
H7 A3970 
[Ne III] A4011 
He I A4026 
[S ii] A4068 
[S ii] A4076 
H(5 A4101 
He I A4120 
He I A4143 
H7 A4340 
[O III] A4363 
He I A4387 
He I A4471 
[Fe III] A4658 
He II A4685 
H^ A4861 
He I A4921 
[O III] A4958 
[O III] A5006 
He I A5015 
NI A5197 
O I A5577 
[N ii] A5754 
He I A5875 
O 1] A6300 
S III] A6312 
O 1] A6363 
N ii] A6548 
Ho A6562 
[N ii] A6583 
He I A6678 
[S ii] A6716 
[S ii] A6730 
He I A7065 
[Ar III] A7135 
He I A7281 
[O ii] A7319 
[O ii] A7330 


-35.9-1-57.7 

0.020±0.001 

0.475±0.009 

0.565±0.011 

0.025±0.001 

0.034±0.002 

0.045±0.001 

0.054±0.001 

0.007±0.001 

0.070±0.001 

0.004 

0.050±0.001 

0.105±0.003 

O.OOSitO.OOl 

0.156±0.004 

O.OlliO.OOl 

0.007±0.001 

0.003±0.001 

0.247±0.005 

0.003±0.001 

0.004±0.001 

0.466±0.009 

0.004 

0.006±0.001 

0.024±0.001 

0.002±0.001 

0.007±0.001 

1.000±0.020 

0.007±0.001 

0.033±0.001 

0.101±0.002 

0.015±0.001 

0.009±0.001 

0.003±0.001 

0.004±0.001 

0.074±0.001 

0.015±0.001 

0.003±0.001 

0.004±0.001 

0.321±0.006 

3.095±0.062 

1.008±0.020 

0.025±0.001 

0.313±0.006 

0.225±0.005 

O.OlOitO.OOl 

0.028±0.001 

0.004±0.001 

0.004 

0.004±0.001 


+49.8+48.7 

0 . 022 + 0.001 

0.694+0.018 

0.959+0.020 

0.027+0.001 

0.037+0.003 

0.049+0.003 

0.058+0.003 

0.016+0.003 

0.082+0.003 

0.012+0.002 

0.052+0.003 

0.111+0.005 

0.013+0.002 

0.009 

0.163+0.007 

0.004 

0.019+0.002 

0.011+0.001 

0.008+0.001 

0.271+0.005 

0.012+0.001 

0.003 

0.464+0.009 

0.008 

0.008+0.001 

0.038+0.001 

0.005+0.001 

0.021+0.001 

1.000+0.020 

0.013+0.001 

0.169+0.003 

0.507+0.010 

0.020+0.001 

0.008+0.001 

0.004 

0.006+0.001 

0.104+0.002 

0.018+0.001 

0.007+0.001 

0.005+0.001 

0.342+0.007 

3.029+0.061 

1.055+0.021 

0.034+0.001 

0.234+0.005 

0.192+0.004 

0.020+0.001 

0.069+0.001 

0.005+0.001 

0.013+0.001 

0.010+0.001 


-73.1-27.3 

0.012+0.001 

0.492+0.010 

0.642+0.013 

0.015+0.001 

0.015+0.001 

0.035+0.001 

0.034+0.001 

0.003+0.001 

0.064+0.001 

0.006+0.002 

0.082+0.002 

0.065+0.002 

0.008+0.001 

0.069+0.008 

0.098+0.003 

0 . 010 + 0.001 

0.011+0.001 

0.006+0.001 

0.275+0.006 

0.004+0.001 

0.481+0.010 
<~ 0.004 
0.007+0.001 
0.021+0.001 
0.006+0.001 
0.010+0.001 
1.000+0.020 
0.008+0.001 
0.049+0.001 
0.144+0.003 
0.013+0.001 
0.009+0.001 
<~ 0.002 
0.001 

0.069+0.001 

0.019+0.001 

0.004+0.001 

0.007+0.001 

0.307+0.006 

3.103+0.062 

0.957+0.019 

0.022+0.001 

0.401+0.008 

0.289+0.006 

0.013+0.001 

0.026+0.001 

0.003+0.001 

0.007+0.001 

0.003+0.001 


-76.2+22.9 

0.018+0.003 

0.507+0.022 

0.658+0.022 

0.022+0.003 

0.027+0.008 

0.036+0.008 

0.048+0.008 

<^^ 0.022 

0.068+0.007 

0.023 

0.046+0.007 
0.093+0.015 
<~ 0.021 
<~ 0.019 
0.138+0.022 

<~ 0.014 
<~ 0.012 
<~ 0.012 
0.253+0.005 
<~ 0.012 
<~ 0.011 
0.465+0.009 
<~ 0.011 
0.010 

0.019+0.003 
<~ 0.009 
<'^ 0.009 
1 . 000 + 0.020 
<~ 0.013 
0.056+0.004 
0.176+0.004 
<~ 0.012 
<~ 0.012 
<~ 0.015 
0.006+0.001 
0.078+0.002 
0.014+0.002 
0.004+0.001 
<~ 0.005 
0.316+0.006 
3.063+0.061 
0.987+0.020 
0.023+0.002 
0.288+0.006 
0.204+0.004 
0.009+0.001 
0.034+0.001 
0.004 

0.007+0.001 

0.004+0.001 


-36.8-73.4 

0.015+0.001 

0.472+0.009 
0.778+0.016 
0.019+0.001 
0.015+0.003 
0.027+0.003 
0.041+0.003 
<~ 0.007 
0.079+0.002 
<~ 0.007 
0.094+0.002 
0.079+0.005 
0.006 
<~ 0.056 
0.116+0.007 
<~ 0.005 
0.011+0.001 
0.004 

0.276+0.006 

0.004+0.001 

0.003 

0.488+0.010 

0.004 

0.011+0.001 

0.018+0.001 

0.003+0.001 

1.000+0.020 

0.008+0.001 

0.051+0.001 

0.155+0.003 

0.013+0.001 

0.015+0.001 

0.003 

0.005+0.001 

0.084+0.002 

0.020+0.001 

0.001 

0.007+0.001 

0.312+0.006 

3.105+0.062 

0.977+0.020 

0.023+0.001 

0.423+0.008 

0.306+0.006 

0.010+0.001 

0.025+0.001 

0.003+0.001 

0.008+0.001 

0.007+0.001 


+68.5+53.4 

0.012+0.002 

0.626+0.013 
0.784+0.016 
0.015+0.003 
0.019+0.006 
0.023+0.006 
0.031+0.005 
0.024+0.005 
0.075+0.005 
0.017+0.003 
0.080+0.005 
0.060+0.010 
0.013+0.004 
<~ 0.076 
0.087+0.015 

<'^ 0.008 
0.014+0.003 
0.008 

0.288+0.006 

<~ 0.006 
0.478+0.010 
<~ 0.006 
0.014+0.002 
0.022+0.002 

0.006+0.002 

1.000+0.020 

0.006+0.002 

0.112+0.002 

0.337+0.007 

0.009+0.002 

0.012+0.002 

0.008+0.001 

0.088+0.002 

0.022+0.001 

0.007+0.001 

0.007+0.001 

0.309+0.006 

3.074+0.061 

1.025+0.020 

0.027+0.001 

0.333+0.007 

0.253+0.005 

0.017+0.001 

0.057+0.001 

0.004+0.001 

0.011+0.001 

0.011+0.001 


+81.6-32.3 

0.015+0.001 

0.352+0.007 

0.580+0.012 

0.019+0.001 

0.017+0.001 

0.036+0.001 

0.041+0.001 

0.010+0.001 

0.067+0.001 

0.003 

0.085+0.002 

0.080+0.002 

0.003 

0.020 

0.120+0.003 

0.008+0.001 

0.017+0.001 

0.002 

0.250+0.005 

0.010+0.001 

0.002+0.001 

0.474+0.009 

0.003 

0.010+0.001 

0.028+0.001 

0.022+0.001 

1.000+0.020 

0.010+0.001 

0.047+0.001 

0.154+0.003 

0.015+0.001 

0.003+0.001 

0.004+0.001 

0.001 

0.101+0.002 
0.008+0.001 
0.003+0.001 
0.003+0.001 
0.294+0.006 
3.182+0.064 
0.898+0.018 
0.031+0.001 
0.243+0.005 
0.170+0.003 
0.014+0.001 
0.036+0.001 
0.004+0.001 
<~ 0.005 
0.004 


[Ar III] A7751 

0.008+0.001 

0.017+0.001 

0.008+0.001 

0.008+0.001 

0.008+0.001 

0.013+0.001 

0.011+0.001 

P18 A8437 

0.004+0.001 

0.004+0.001 

0.004+0.001 

0.019 

0.006+0.001 

0.011+0.002 

0.007+0.001 

O I A8446 

0.006+0.001 

0.007+0.001 

0.011+0.001 

0.020 

0.007+0.001 



P17 A8467 

0.004+0.001 

0.005+0.001 

0.008+0.001 

<~ 0.019 

0.004+0.001 

<~ 0.005 

0.006+0.001 

P16 A8502 

0.006+0.001 

0.006+0.001 

0.009+0.001 

0.020 

0.005+0.001 

<~ 0.005 

0.005+0.001 

P15 A8545 

0.007+0.001 

0.007+0.001 

0.012+0.001 

<~ 0.019 

0.006+0.001 

<~ 0.005 

0.007+0.001 

P14 A8598 

0.007+0.001 

0.008+0.001 

0.013+0.001 

0.020 

0.007+0.001 

0.008+0.002 

0.008+0.001 

P13 A8665 

0.009+0.001 

0.011+0.001 

0.017+0.001 

<~ 0.021 

0.009+0.001 

0.008+0.002 

0.010+0.001 

N I A8683 

0.002+0.001 

0.004 

0.006+0.001 

<~ 0.021 

0.002 

<~ 0.005 

0.002 

P12 A8750 

0.005+0.001 

0.013+0.001 

0.012+0.001 

<~ 0.021 

0.020+0.001 

0.037+0.002 

0.012+0.001 

Pll A8862 

0.010+0.001 

0.019+0.001 

0.021+0.001 

0.022 

0.017+0.001 

0.028+0.002 

0.023+0.001 

PIO A9015 

0.024+0.001 

0.022+0.001 

0.024+0.001 

0.021+0.002 

0.024+0.001 

0.017+0.001 

0.025+0.001 

[S III] A9068 

0.238+0.005 

0.367+0.007 

0.192+0.004 

0.245+0.005 

0.167+0.003 

0.293+0.006 

0.284+0.006 

P9 A9229 

0.030+0.001 

0.030+0.001 

0.032+0.001 

0.031+0.003 

0.025+0.001 

0.029+0.001 

0.036+0.001 

[S III] A9530 

0.644+0.013 

0.969+0.019 

0.554+0.011 

0.699+0.015 

0.479+0.010 

0.837+0.017 

0.788+0.016 

P8 A9546 

0.020+0.001 

0.044+0.001 

0.025+0.001 

0.019+0.005 

0.054+0.001 

0.030+0.001 

0.061+0.001 

Gif/3 

0.297+0.013 

0.450+0.013 

0.221+0.013 

0.250+0.013 

0.433+0.013 

0.542+0.013 

0.148+0.013 

Ffr/3 

257.9+0.4 

715.5+2.9 

391.8+0.3 

202.3+1.2 

190.2+0.3 

77.9+0.2 

297.6+0.2 

EW HP 

202.3 

139.8 

102.6 

143.3 

121.6 

90.9 

114.5 

EWif^ 

1331.5 

1109.7 

647.0 

826.5 

910.7 

666.0 

863.8 


Note. — Emission line fluxes. Detections at a significance of less than 3 a are given as upper limits. Full table available 
online only. 
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interactions such as relativistic corrections. However, 
they note that with the exception of [S II], the result¬ 
ing theoretical temperature dispersions have been found 
to be no larger than 10%. Followin g the recommen¬ 
dations of IMendoza fc Bautista] (120141 . and considering 
the discussion o f effect ive collision strengths for [O III] 
by iStorev et al.l (|20l41 . we adopt the most up-to-date 
atomic data (see Table |4]). 

To demonstrate the significance of atomic data up¬ 
dates, in Figure [3] we present a comparison of temper¬ 
atures determined for the K03 MlOl r egions using (1) 
the at omic data as in K03: specihca lly, iTaval fc Guotal 
(|1999ll and iKauffman fc Sugad (1198611 for [S 111], and (2) 
the current database adopted here (see Table Q. Tem¬ 
peratures determined with the old atomic data are plot¬ 
ted as open circles, while the updated values are closed 
circles, and the set of points belonging to the same H 11 
region are connected. For reference, the one-to-one re¬ 
lation is plotted as a solid line and the theoretical re¬ 
lationship of G92 is shown as a short dashed line. The 
comparison in Figure [3] shows a shift of roughly 10% in 
T[S 111] (^1000 K) in the sense that using the newer 
atomic data results in lower values of T[S III]. Addition¬ 
ally, in Figure U we plot the difference in electron tem¬ 
perature due to updating the atomic data versus the new 
electron temperature. The old temperatures are deter¬ 
mined using the same atomic data as K03, whereas the 
new temperatures use the updated atomic data listed in 
TableHl The temperature trends shown in Figure|3]show 
that the atomic data for [N II] and [O III] have remained 
fairly stable, whereas the data for [O II] and [S III] have 
not yet converged. As we will see in the next section, 
our observations show significant discrepancies between 
values of T[S III] and T[0 III], in the sense that some 
T[0 III] values are signihcantly higher. 

3.3. Electron Temperature and Density Determinations 

To establish a low-uncertainty electron temperature es¬ 
timate, we define a quality criterion requiring a 3 (t de¬ 
tection of the auroral line measurement (see § 12.31) . The 
high throughput of the MODS over a large optical and 
near-UV wavelength range allows us to take advantage 
of the following auroral lines: [S II] AA4068,4076, [O III] 
A4363, [Ar III] A5192, [N II] A5755, [S III] A6312, and 
[O II] AA7320,7330. In the NGC 628 spectra, there are 
very few [S II] AA4068,4076 and [Ar III] A5192 detec¬ 
tions, so we will focus on the four lines with the most 
numerous detections. We detect the following auroral 
line measurements at a strength of 3cr or greater: [O III] 
A4363 in 18 H II regions, [N II] A5755 in 29 H II regions, 
[S III] A6312 in 40 H II regions, and [O II] AA7320,7330 
in 40 H II regions. Due to the previo us concerns with 
[O II] temperature measurements ('e.g.. IKennicutt et al.l 
I2003bll . we do not report abundances calculated from 
this diagnostic. Rather, we defer to a future analysis of 
the entire GHAOS sample. The result is a sample of 45 
H II regions with direct abundances from [O III], [N II], 
or [S III]; where 40 of these regions have multiple tem¬ 
perature determinations. The auroral lines measured for 
each H II region are listed in Table [3] 

The [S II] AA67I7,6731 ratio is used to determine the 
electron densities, which are conhrmed with density de¬ 
terminations from the [O II] AA3726,3728 line ratio 
For all densities which are consistent with the low den- 



T[0 III] (xi0'‘ K) 

Fig. 3.— The effect of updating atomic data as demonstrated 
by comparisons of [O III] and [S III] temperature determinations 
for the H II regions in MlOl from K03. Te mperatures determined 
wit h the atomic dataset of KQ3 ('specificallv. lTaval &; Guptal II1999I) 
and IKauhman &: Su^ B IToM) for [S III]) are plotted as open cir¬ 
cles, whereas values determined using the sources adopted here (see 
Table [JJl are closed circles. Sets of points belonging to the same 
H II region are connected. For reference, the equivalence relation 
is plotted as a solid line and the theoretical relationship of G92 is 
shown as a short dashed line. 

sity limit, abundance calculations assume Ue = 10^ cm“^ 
(which is consistent with the la upper bounds and pro¬ 
duces identical results for all lower values of rie), else the 
derived higher electron density (ue > 10^ cm“^) is used. 
The electron temperature and density determinations are 
listed in Table [5] (full table is available online). 

3.4. Comparing Electron Temperature Measurements 

Using the 45 H II regions in our sample with multi¬ 
ple auroral line measurements, all of which were dou¬ 
ble checked by hand, we compare the calculated electron 
temperatures in Figure [SJ The red dashed lines show the 
line of equality and the black short dashed lines are the 
expected relationships from the photoionization models 
calculated by G92. For each set of variables we calculate 
the scatter in the data, with errors in both coordinates, 
around the line of equality and determine how much can 
be attributed to a physi cally significant intrin sic scatter 
using the MPFITEXY (|Williams et al.llMlFI routine 
in IDL. The intrinsic scatter is determined by minimiz¬ 
ing the departure of (degrees of freedom) from 1 by 
introducing an additional scatter term to the weighting 
of each data point. Additionally, in two cases we also 
perform a least-squares fit, which is shown with a solid 
line. The calculated total and intrinsic scatters, at and 

At the resolution of MODS, the [O II] AA3726, 3728 doublet is 
partially blended for all observations, and so measuring this ratio 
serves primarily as a sanity check of the [S II] density determina¬ 
tion. 
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TABLE 4 

Sources of Updated Atomic Data for CHAOS 


Ion 


Radiative Transition Probabilities 


Collision Strengths 



j^isieliu^^^aL ^009) 

^tore^^^a^ 

Taval (20U) 

^cLaughlii^^^ell (200^ 
^unoz^^urgo^^e^^^ ^^0^ 
^^a^^^atsarinw 

^^Hudsoir^^aI^2Q^^"^ 


Note. — The atomic data used by the 5-level atom approximation has been 
updated with current values from t he lite r ature as reported here. 

* Agrees with updated values from [Tavall i2011f) . 

^ Equivalent to [ Tachiev &c Froese Fische^ II2002I) . as recommended by 
IStasihska et all ||201% . 


TABLE 5 

Ionic and Total Abundances for MODS/LBT Observations of NGC 628 


Ha; Region 

-35.9-1-57.7 

-1- 49.8-148.7 

-73.1-27.3 

-76.2-122.9 

-36.8-73.4 

±68.5±53.4 



9000±100 

7300±100 

7500±500 

8200T100 

9300T200 

^i^S^measured {^) 







t(N2)77ieasiired (K) 

6800±100 

7200±200 


7300±400 

7300T200 

8100T300 

t(S3)measti,reci (^) 

5700±100 

6400±100 

6400±100 

6000±400 


6600T200 

'^e.^measured) (cm 

57 

248 

59 

37 

64 

126 

H02)used (K) 

6800±500 

7200±200 

6900±500 

7300±400 

7300T200 

8100T300 

tiOsiused (K) 

4800±500 

5700±500 

5600±500 

5200±500 

6100T500 

6000T500 

t(N2)used (K) 

6800±200 

7200±200 

6900±500 

7300±400 

7300T200 

8100T300 

t(S3)used (K) 

5700±200 

6400±200 

6400±200 

6000±400 

6800T500 

6600T200 

rieiused) (cm-6) 

100 

248 

100 

100 

100 

126 

0-k/H+ (106) 

28.1±1.4 

32.5±2.0 

26.9±5.3 

20.4±3.2 

21.4±1.3 

14.2±1.2 

0-k+/H-k (105) 

11.7±3.8 

22.1±5.2 

6.58±1.55 

12.9±3.9 

4.44±0.90 

11.3T2.4 

12 + log(0/H) (dex) 

8.60±0.04 

8.74±0.04 

8.52±0.07 

8.52±0.07 

8.41±0.03 

8.41±0.05 

N-I-/H+ (106) 

64.8±2.1 

56.5±2.2 

57.2±6.9 

49.1±4.8 

49.2±1.8 

35.2±1.8 

N ICE 

1.417±0.161 

1.680±0.200 

1.245±0.318 

1.635±0.361 

1.207T0.101 

1.796T0.240 

log(N/0) (dex) 

-0.64±0.03 

-0.76±0.03 

-0.67±0.10 

-0.62±0.08 

-0.64±0.03 

-0.61T0.04 

12 -1 log(N/H) (dex) 

7.96±0.05 

7.98±0.05 

7.85±0.12 

7.90±0.10 

7.77T0.04 

7.80T0.06 

S-k/H+ (lO'O 

43.1±1.4 

29.2±1.1 

51.4±6.0 

30.7±2.9 

46.1T1.6 

26.7T1.3 

S-k+/H-k (IQl') 

178.±5.4 

188.±4.4 

104.±2.2 

160.±17. 

75.6T7.4 

141.±5.3 

S ICE 

1.083±0.108 

1.151±0.115 

1.026±0.103 

1.142±0.114 

1.013±0.101 

1.172±0.117 

log(S/0) (dex) 

-1.22±0.06 

-1.34±0.06 

-1.32±0.09 

-1.19T0.09 

-1.32±0.06 

-1.11±0.06 

12 + log(S/H) (dex) 

7.38±0.05 

7.40±0.04 

7.20±0.05 

7.34±0.06 

7.09T0.05 

7.29±0.05 

Ne+-k/H+ (106) 


24.8±0.58 

14.1±2.8 



26.9±0.63 

Ne ICE 

3.398±1.144 

2.471±0.628 

5.085±1.461 

2.575T0.878 

5.823T1.224 

2.256±0.535 

log(Ne/0) (dex) 


-0.95±0.10 

-0.67±0.13 



-0.62±0.09 

12 -1- log(Ne/H) (dex) 


7.79±0.11 

7.85±0.15 



7.78±0.10 

Ar-l+/H-k (10') 

14.7±0.53 

21.5±0.56 

9.08±0.22 

13.3T2.0 

7.02±0.83 

14.9±0.69 

Ar ICE 

1.916±0.192 

1.569±0.157 

2.730±0.273 

1.602±0.160 

3.140T0.314 

1.510±0.151 

log(Ar/0) (dex) 

-2.15±0.06 

-2.21±0.06 

-2.13±0.08 

-2.19T0.10 

-2.07T0.07 

-2.05±0.07 

12 + log(Ar/H) (dex) 

6.45±0.05 

6.53±0.05 

6.39±0.05 

6.33T0.08 

6.34±0.07 

6.35±0.05 


Note. — Electron temperatures and ionic and total abundances for objects with an [O III] A4363, [N ll] A5755, 
or [S III] A6312 line signal to noise ratio of 3cr or greater. Electron temperatures were calculated using the [O III] 
(A4959 -h A5007)/A4363, [N ii] (A6548 -h A6584)/A5755, or the [S ni] (A9069 -h A9532)/A6312 diagnostic line ratio. Full 
table available online only. 


(7i respectively, are displayed in Figure [5j 

The left column of Figure[5]shows T[N II], T[0 III], and 
T[0 II] as a function of T[S III], The top panel of this 
column shows a very tight correlation between T[N II] 
and T[S III]. The dispersion in this relationship about 
the line of equality is 240 K, and is consistent with be¬ 
ing due entirely to observational uncertainties (i.e., very 
little internal scatter). The middle panel shows a large 
scatter for the comparison of T[0 III] and T[S III], with a 

16 MPFITEXY is dependent on the MPFIT package IIMarkwardtl 

[2003). 


dispersion around the line of equality of 1270 K, indicat¬ 
ing an intrinsic dispersion of 1160K. The bottom panel 
shows a comparison of the T[0 II] and T[S III] measure¬ 
ments, where the scatter is, again, quite large at 1610 K, 
indicating an intrinsic dispersion of 1560 K. This panel 
is similar to the T[0 III] versus T[S III] comparison in 
having both a large dispersion and a systematic offset to 
smaller T[S III] values. 

The right column of Figure [5] compares T[N II], 
T[0 III], and TJO II] to one another. The top panel 
shows T[N II] versus T[0 II]. In this comparison, the 










































11 


o 

X 


o 

I 

cu 

c 


“T-1- 1 -1-1-1- 1 - 1 - 1 -1-1- 1 -!— 








T[oin] 

T[NII] 

T[OII] 

T[SIII] 


1 t t 1 1 I 


t 1 1 1 t 


0.6 0.8 1 1.2 
T(new) (xlO^ K) 


1.4 


Fig. 4.— The effect of updating atomic data as demonstrated 
by the difference of new and old electron temperatures versus new 
temperature. Data is color coded for each ion: T[0 III] - red, 
T[N II] - orange, [O II] - green, and T[S III] - blue. Little difference 
is seen for T[N II] and T[0 III], whereas significant updates have 
been made to the [S III] and [O II] atomic data and resulting 
temperatures. 


scatter is again large, with a dispersion of 1550 K, in¬ 
dicating that essentially all of the dispersion is intrin¬ 
sic (formally 1510 K). In this comparison, where roughly 
equal temperatures are expected from photoionization 
models, the data appear to fall into two groups. Roughly 
half of the data are shifted from the equality relationship 
by ~ 1400 K toward higher T[0 11]. The other half of the 
data show an even larger offset in the sense that the 
T[0 11] are hotter than the T[N II] estimates. Of the 
outliers, three fall into a group where the average off¬ 
set is close to 4000 K. Most of the very large dispersion 
is due to this triad of outlying points. Overall, T[0 II 
is systematically larger than both T[N II] and T[S III 
as noted in previou s studies (e.g., iPilvuein et al.ll2005 
lEsteban et al.ll^OOQf) . The middle panel shows T[0 III] 
versus T[N II]. Although there are only eight points in 
this comparison, all but one of the measurements are 
in agreement with the expected relationship from G92. 
Finally, in the bottom panel, the comparison between 
T[0 II] and T[0 III] shows a large dispersion of 1140 K, 
of which 1020 K is due to intrinsic scatter. These points 
show a similar trend to that seen in the T[0 11]—T[S 111] 
comparison, with preferential scatter to higher values of 
T[0 11]. 

The magnitude of the scatter in the T[0 III]—T[S III] 
relationship came as a surprise. In fact, the CHAOS 
program was designed around the assumption that a se¬ 
cure measurement of T[0 III] implied a secure abundance 
(with the understanding that there could be small bi¬ 
ases d ue t o temperature inhomogeneit ies). iHaeele et al.l 
(|2006f) and iPerez-Montero et al.l (120061) have both shown 
a poor correlation between T[0 III] and T[S III], and 
iBinette et al.l (|2012D showed definitively that discrepan¬ 


cies between T[0 III] and T[S III] are very common in the 
sense that T[0 III] is higher than T[S III]. We showed in 
Section l3^ that recent updates to atomic data enhances 
the T[0 III]-T[S III] discrepancy by uniformly shifting 
[S III] temperatures cooler, but cannot be entirely re¬ 
sponsible for the magnitude of discrepancies observed. 
Although many of the points in NGC 628 are consistent 
with the equality trend or the relationship derived by 
G92, there are a significant n umber of outli e rs, in agree¬ 
ment with the trend noted bv IBinette et al.l (|2012h . This 
large number of outliers, and the similar trends noted in 
the literature, call into question the reliability of [O III] 
as a temperature indicator in H II regions. We discuss 
the discrepant temperature measurements further in Sec¬ 
tion [5T] 

In sum, the T[N II] and T[S III] measurements appear 
to be generally consistent, especially for Tg < lO"* K, 
while the T[0 III] and T[0 II] measurements are often 
discrepant. Because a significant fraction of the T[0 III] 
and T[0 II] measurements appear to be biased to higher 
than expected values, abundances based on the these 
temperature measurements exhibit both larger disper¬ 
sions and systematic offsets. These results guide our cho- 
sen methodology for deriving abundances in Section 14.11 

4. A HOMOGENEOUS ABUNDANCE ANALYSIS BASED ON 
[S III] AND [N II] TEMPERATURES 

4.1. Methodology for Ionic and Absolute Abundance 
Calculations 

From the comparisons in Section 13.41 it is clear that 
there are significant temperature measurement discrep¬ 
ancies. Because T[N II] and T[S III] show a very good 
correlation at low temperatures, we interpret this to 
mean that these temperatures are free from whatever sys¬ 
tematic effects are affecting the T[0 III] and T[0 II] mea¬ 
surements. Fortunately, the vast majority of our sample 
have measurements of T[N II] or T[S III] (or both). Thus, 
we calculate absolute and relative abundances based pri¬ 
marily on temperatures determined from T[N II] and/or 
T[S III] in conjunction with the G92 scaling relation¬ 
ships for the other ionization zones. Specifically, when 
T[N II] and T[S III] are both measured, T[N II] is used 
for the low ionization zone, T[S III] is used for the inter¬ 
mediate zone, and T[S III] in combination with the G92 
relationship is used for the high ionization zone. In the 
cases with only one temperature measurement, the tem¬ 
peratures in the other zones are inferred from the G92 
relationships. Since we have [S III] or [N II] temperatures 
for 43 of the 45 H II regions, we can adopt this relatively 
homogeneous method for nearly the entire dataset. We 
base the subsequent analysis on this prioritization. For 
completeness, we present an abundance analysis based 
on T[0 III] in Appendix lAl 

Ionic abundances relative to hydrogen are calculated 
using: 

N{x^) r.. 

N{H+) Ihp JAP) ■ ^ ^ 

The emissivity coefficients, jA(i)j which are functions of 
both temperature and density, are determined using a 5- 
level atom approximation with the updated the atomic 
data reported in Table [U 

Total oxygen abundances (0/H) are calculated from 
the simple sum of and The other 
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Fig. 5. — Comparisons of the relations between the [O III], [N II], [O II], and [S III] Tg temperatures for H II regions with more than one 
temperature measurement. The red long-dashed lines represent equal temperatures and black short-dashed lines are the G92 relationships. 
For the T[N II] versus T[S III] and T[0 III] versus T[S III] panels, the error-weighted least-squares best fit is shown by solid lines. The total 
and intrinsic scatter around the line of equality is give for each comparison. 
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abundance determinations require ionization correction 
factors (ICF) to account for unobserved ionic species. For 
nitrogen, w e employ the co mmon assumption that N/0 
= N+/0+ (jPeimbertIIl967li . This allows us to directly 
compare our resu lts with other studies in the literature. 

I Nava etldl ()2006[1 have investigated the validity of this 
assumption, and concluded that it is valid at a precision 
of about 10%. 

Collisionally excited emission lines of sulfur, neon, and 
argon are also observed in many of our spectra. For 
Ne we use a fai rly straightforward IC F: ICF(Ne) = (0+ 
+ 0++)/0++ (j Crockett et al. 1 1^00611 . S and Ar present 
more complicated situations as both and S’*"'*'’*' lie in 
the 0 +~*~ zone and Ar~*~+ sp ans both the 0+ and the 0++ 
zones. iThuan et al.l (|1995r) provide analytic ICF approx¬ 
imations for both S and Ar using the model calcul ations 
of photoionized H II reg ions bv iStasinskal (119901) . We 
employ these ICFs from iThuan et al.l ( I995f) to correct 
for the unobserved S+^, Ar+^, and Ar+'^ states: 

ICF(S) = , ^ , 

^ ' S+-HS++ 

= [0.013 + x{5.10 + a;[-12.78 + x(14.77 - O.lla:)]}] 

(4) 

Ar 

ICF(Ar) = 

= {0.99 -I- 3;[0.091 -I- a;(-1.14 + 0.077a:)]}-l (5) 

Ar 

= —XX = [0'15-|-a;(2.39-2.64a:)]"\ (6) 

where x = QA/O and the second Ar equation is used 
when [Ar IV] A4740 emission is not observed. Ionic 
and total abundances are listed for 45 H II regions in 
NGC 628 in Table El 

4.2. Alpha-Element Relative Abundances 

a-elements are thought to be produced in massive stars 
such that sulfur, neo n, and argon are produce d in paral¬ 
lel with oxygen fe.g.. lWooslev fc Weaverl[l995ll . Figure El 
shows the sulfur, neon, and argon abundances measured 
for NGC 628 versus oxygen abundance. Error-weighted 
least-squares fits produce gradients within one standard 
deviation of a flat (zero) slope, and so we adopt the 
error-weighted average alpha-element abundances as con¬ 
stants, which are shown in FigureElas shaded boxes span¬ 
ning one standard deviation. The sample covers a range 
of 8.0 < 12-|-log(0/H) < 9.0 with no obvious trends in 
sulfur, neon, or argon abundance. S/0, Ne/0, and Ar/0 
all show very small intrinsic dispersions of 0.08, 0.08, 
and 0.02 dex, respectively. We report average a-element 
to oxygen abundance ratios in Table El In general, we 
find that the relative abundances of the a-elements in 
NGC 628 behave as expected for the nucleosynthetic 
products of massive stars. The constant values over the 
range in metallicity are consistent with a universal value 
for the upper mass IMF. This is in agreem ent with the 
recent study of the IMF in MSI clusters bv lWeisz et al.l 

(IM^ . 

4.3. Radial Gradients in N/0 and 0/H 

We display the derived N/0 and 0/H abundances in 
Figure [71 as a function of galactocentric radius. Because 
the locations of individual H II regions are known with 



Fig. 6. — a-elements are plotted versus oxygen abundance for 
NGC 628. The error-weighted average a-element abundance is 
illustrated in each panel by a shaded box centered on the average 
and extended to dilcj. The measured dispersions and calculated 
intrinsic dispersions are indicated in the upper right of each panel. 
All three panels demonstrate constant a-element/O ratio trends 
with small dispersions. 

high precision relative to one another, we consider only 
the uncertainties associated with oxygen abundance 
here. We characterize the N/0 observations with a 
single, linear, error-weighted least-squares fit (solid line): 

log(N/0) 

= (-0.523 ± 0.034) + (-0.077 ± 0.006) x R (dex/kpc) 

= (-0.521 ± 0.035) -I- (-0.849 ± 0.064) x R (dex/fl25), (7) 

with a measured dispersion in log(N/O) of ct = 0.083 dex. 
Although this dispersion is small, the statistical analysis 
shows that most of this dispersion is intrinsic (di = 0.060 
dex). 

The bottom panel of Figure [7] shows the values of 
0/H as a function of galactocentric radius. The re¬ 
sulting best fit to the 45 objects in the sample is given by: 

12 -k log(0/H) 

= (8.834 ± 0.069) -|- (-0.044 ± 0.011) x R (dex/kpc) 

= (8.835 ± 0.069) -|- (-0.485 ± 0.122) x R {dex/R 25 ), (8) 

with a dispersion in log(0/H) of cr = 0.165 dex. Ac¬ 
counting for the observational uncertainties results in an 
intrinsic dispersion of 0.159 dex. This is significantly 
larger than the dispersions seen in the relative abun¬ 
dances. Since it is conceptually difficult to produce sig¬ 
nificant dispersions in absolute abundances without sig¬ 
nificant dispersions in relative abundances (due to time 
delays of delivery and mixing timescales), it is possible 
that this relatively large dispersion could be due to dis¬ 
crepant temperatures. The final CHAOS gradients are 
tabulated in Table El 
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TABLE 6 

Adopted Abundance Gradients for NGC 628 


y 


X 


Equation of Correlation 


<^i 

12+log(0/H) 

VS. 

R (dex/kpc) 

y = 

(8.834±0.069) -1- (-0.044±0.011)xX 

0.165 

0.159 



R (dex/it25) 

y = 

(8.835±0.069) -|- (-0.485±0.122)xX 



log(N/0) 

vs. 

R (dex/kpc) 

y = 

(-0.523±0.034) -1- (-0.077±0.005)xX 

0.083 

0.060 



R {dex/i?25) 

y = 

(-0.521±0.035) -1- (-0.849±0.064)xX 



log(S/0) 

vs. 

12-Hlog(0/H) 

y = 

-1.413 

0.108 

0.082 

log(Ne/0) 

vs. 

12-Hlog(0/H) 

y = 

-0.714 

0.125 

0.091 

log(Ar/0) 

vs. 

12-Hlog(0/H) 

y = 

-2.219 

0.082 

0.038 


Note. — The adopted best fits to the abundance gradients measured for the 3 cr NGC 628 
CHAOS data. 



Fig. 7.— N/O and O/H plotted versus galactocentric radius for 
NGC 628. The solid lines represent the least-squares linear fits to 
the data when considering the errors in N/O and O/H abundance 
respectively. 

5. DISCUSSION 

5.1. Discrepant Temperature Measurements 
5.1.1. Suspect [O III] Temperatures 

Historically, [O III] A4363 is the best studied direct 
temperature auroral line, owing to its relative ease of 
observation, abundance of emitting ions, and increas¬ 
ing strength in the moderate- to low-metallicity H 11 
regions of nearby dwarf and (all but the central re¬ 
gions of) spiral galaxies. Other common temperature- 
sensitive ratios, such as [N II] I(AA6548,6584)/I(A5755), 
[O II] I(AA7320,7330)/I(AA3726,3628), and [S III] 
I(AA9069,9532)/I(A63I2), are limited by their weaker 
abundances, origination in the smaller low-ionization 
zone, and combination of optical and near-infrared line 
ratios. Because of these factors, the intensity ratio [O III] 
A4363/AA(4959-|-5007) has long been the chief diagnostic 
of nebular gas temperature. 

In theory, the temperature measurements shown in 
Figure El should fall along the lines of the relationships 
from the photoioniz ation models of G92, but in many re¬ 
gards they do not. iKennicutt et al.l (|2003bD raised con¬ 


cerns about T[0 II], and additional studies have pointed 
out other significant temperature discrepancies (e.g., 
Breso in et al.l 120051 iPilviigin et al1l2009t lEsteban et alJ 

20091 : Ziirita fc BresolinI 120121 : iBinette et al.ll2012t) . For 
the CHAOS data. Figure O demonstrates a surprisingly 
large dispersion for T[S III] versus T[0 III], yet a tight 
relationship for T[N II] versus T[S III]. We investigate 
these two relationships further here. 

We reproduce the CHAOS T]N II]— T[S III] data in 
the left panel of Figure [H As mentioned earlier, the 
dispersion in this relationship about the line of equal¬ 
ity is small. K03 not did see such a tight relationship, 
but there were only 7 H II regions with measurements 
of both T[N II] and T[S III]. The K03 observations have 
been added to Figure [5] for comoarisorFI. For consis¬ 
tency, the K03 temperatures are recalculated using the 
K03 emission line fluxes and our adopted atomic data. 
The least-squares fit to the combined set sees an in¬ 
crease in the dispersion due to the additional scatter at 
higher temperatures. A very tight relationship between 
these two temperature diagno s tics has been noted be¬ 
fore (e.g.. iBresolin et al.l 120051 : IZurita k, Bresoliiil 120121 : 
IBinette et al.ll2012l l. The emergence of an excellent cor¬ 
relation indicates that T[N II] and T[S III] provide highly 
consistent measurements of H II region temperatures and 
increases confidence in their reliability at low tempera¬ 
tures. 

In the right panel of Figure |8] we plot the CHAOS 
T[0 III]—T[S III] measurements and compare them to 
the K03 data. K03 found excellent agreement between 
T[0 III]— T[S III] measurements and the expectation 
from the photoionization modeling of G92. Our observa¬ 
tions of NGC 628, coupled with updated atomic data, are 
at odds with the MIDI study, showing a significant dis¬ 
persion in the T[0 III]—T[S III] comparison. The least- 
squares ht to the combined data set results in a total 
dispersion of 1090 K which corresponds to an intrinsic 
uncertainty of 950 K. As shown in Figure |3l the original 
K03 (S III] temperatures are higher by about 10%. The 
result is a shift in the T[0 III]— T(S III] relationship away 
from the G92 relationship, and is equivalent to the off¬ 
set ob served for the CHAOS data. iPerez-Montero et al.l 
(j2006l l suggest a theoretical T(0 III]— T(S III] relation¬ 
ship that also differs from the G92 standard, identifying 
new atomic coefficients as the source of the difference. 
However, the new atomic data have minimal effect on 
the observed dispersions. One very important thing to 
note from Figure [8] is the difference in the temperature 
ranges covered by the observations in the two plots. The 
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left hand panel concentrates on a tight correlation for 
values of T[N II] below 9,000 K, while the bulk of the 
T[0 III] values in the right hand panel, which show signif¬ 
icant scatter, are above ^ 7,000 K. Thus, the conclusions 
from K03 and the present study could be understood as 
an increasing dispersion in the T[0 III]— T[S III] relation¬ 
ship with decreasing temperature. This result appears to 
be supported by the right hand panel of Figure [H 


5.1.2. Sources of Temperature Discrepancies 

We now turn to the possible causes of the signifi¬ 
cant temperature discrepancies demonstrated by Fig¬ 
ures [5] and [H Due to the moderately-high resolution 
of MODS, [O III] A4363 is resolved from the Hg A4358 
sky line. Additionally, we examine the [O III] auroral 
line emission by hand in both the one-dimensional and 
two-dimensional spectra, eliminating concerns of cosmic 
rays, Bal mer absorptio n , or o ther measurement errors. 
Recently, IBinette et al.l (120121) have highlighted the dis¬ 
crepancy between the T[0 III] and T[S III] tempera¬ 
tures and investigated possible sources of this disagree¬ 
ment: temperature inhomogeneities, metallicity inho¬ 
mogeneities, a non-Maxwell-Boltzmann distribution of 
electron energies, and shocks. Temperature inhomo¬ 
geneities within ind i vidua l H II regions were first pro¬ 
posed by iPeimbe^ (119671) , where non-uniform temper¬ 
ature distributions tend to result in average observed 
temperatures bia sed higher than tru e local tempera¬ 
tures. However, IBinette et al.l (|2012h found that sim¬ 
ple temperature fluctuations are insufficient to account 
for the magnitude of the T[0 III]—T[S III] disparity, 
especially at the low electron temperatures character¬ 
istic of high abundance H II regions (see their Figure 
4). Instead, Binette et al. conclude that meta llicity in¬ 
homogeneities (e.g., iTsamis et al.l 12003 . I20II1) . adopt- 
i ng a K-distribution for the electron energy distribution 
(jNicholls et al.l 120121 ). and pollution by shock waves all 
successfully reproduce the observ ed T[0 III] excesses 
see th eir Figures 5 — 8). However. iMendoza fc Bautistal 
20I4D found K-distributions only account for observed 
temperature variations in low-excitation regions, and do 
not reconcile the differences in high-excitation regions. 

For the low elect ron temperat u res as sociated with 
the CHAOS data, IBinette et al.l (j20I2[) demonstrate 
that plane-parallel shock models which propagate 
within a photoionized gas layer successfully bridge the 
T[0 III]—T[S III] gap. The shocked nebula models show 
a diminishing excess of T[0 III] with decreasing nebular 
metallicity (increasing T) similar to what we see in our 
data, where the largest excesses occur with decreasing 
temperature (increasing metallicity) towards the inner 
disk. By combining the continuous starburst model of 
Kewlev et alJ (l200ll) and the sho c k-|-pre cursor model of 
Allen et al.l ( 20081) . IBinette et al.l (I20I2D suggest an ad¬ 
ditional region which is a mixing zone of the photo- and 
shock-ionization zones. 

In Figure [9] we reproduce the s hock model data 
from Figure 8 of IBinette et al.l ()20I2[) and compare the 
CHAOS data for NGC 628. The solid points are the 
H II regions which have both measured [O III] and [S III] 


In order to assure the quality of the data compared and uni¬ 
form methodology, a 3(t S/N cut in flux was made for the MlOl 
data. 


temperatures, and are color coded by abundance. Points 
in this plot show that (with the exception of one ma¬ 
jor outlier: NGC 628-168.2-I-I50.8) T[0 III] is in excess 
of T[S III], where increased temperatures subsequently 
lead to artificially low oxygen abundances. This idea is 
congruent with the picture shown in Figure [9j the colors 
of the points indicate that their abundances are higher 
than those of the models to which they are closest. If 
their true photoionized T[0 III] is lower (~T[S III] ), 
then the true oxygen abundance is higher and in better 
agreement with the model abundances. 

It is clear from studies in the literature and from our 
CHAOS observations of NGC 628 that emission line spec¬ 
tra for some H II regions are affected by one or more 
processes which lead to discrepant temperatures derived 
from auroral line observations. By comparing the present 
observations with theoretical models, we have shown that 
shocks, as an example, can contribute to discrepant tem¬ 
peratures; however, the true nature of the perturbations 
is not clear. It is interesting to note that the auroral lines 
with the highest excitation energies (e.g., [O III] A4363 
and [O II] AA7320, 7330, with excitation energies of 5.3 
and 5.0 eV respectively) show the largest discrepancies, 
while those with lower excitation energies (e.g., [N II] 
A5755 and [S III]A63I2, with excitation energies of 4.1 
and 3.3 eV respectively) are less sensitive to these per¬ 
turbations. The higher excitation energies of the oxygen 
temperature diagnostic lines may contribute to making 
these lines more vulnerable to the perturbations, regard¬ 
less of the physical cause. This leads quite naturally to 
the result that the biases are larger in higher abundance 
(lower temperature) H II regions, in agreement with the 
observations. From the present NGG 628 dataset, we 
cannot deduce a single likely cause for the temperature 
discrepancies. However, in Appendix IbI we suggest rela¬ 
tive Ar abundances as excellent diagnostics for identify¬ 
ing which spectra are most affected. We will revisit this 
problem with the entire GHAOS dataset. 

5.2. Comparison of Gradients with Measurements from 
the Literature 

Although NGC628 is one of the best-studied spiral 
galaxies, there are only a handful of previous studies 
with direct abundance measurements. We first perform 
a straightforward comparison to these studies, and then 
extend t his comparison to the integra l field spectroscopy 
study bv iRosales-Ortega et al.l (1201 ID . 

Three previous studies present dire c t abu ndance mea¬ 
surements for NGG 628: iBerg et al.l ()20I3L 14 regions, 
hereafter BI31. iGr oxall et al.l (1201317 regions, hereafter 
GI3), and iGastellanos et al.l (j2002L I region), but only 
the first reports nitrogen or a-element abundance ra¬ 
tios. Note that the three studies derive electron temper¬ 
atures using different methods. BI3 used the [O III] and 
[N II] auroral lines to determine electron temperatures, 
iGroxall et al.l (I20I3D used the infrared [O I II] 88^m fine- 
structure line, and [Castellanos et al.l (j2002D found agree¬ 
ment amongst four temperature determinations from the 
[O III], [S III], [N II], and [S II] auroral lines observed in 
one H II region. We have not attempted to make any cor¬ 
rections to place these different abundance derivations on 
a common scale. 

5.2.1. Radial Gradients in N/0 and 0/H 
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Fig. 8. — Comparisons of the relations between T[N II] and T[0 III] v ersus T[S III]. The lines a nd solid points are the same as in Figure!^ 
A parallel comparison is provided by the MlOl direct measurements of IKennicutt et al.l <|20Q3bl'l . where the open blue points correspond to 
the temperature values presented in their work and the closed blue points represent electron temperatures recalculated from their emission 
line fluxes with the adopted updated atomic data. 



Fig. 9.— T[0 III] versus T[S III] plot for NGC 628 in comparison 
to models including the effects of shocks from IBinette et al.l <I2012IV 
Abundances are divided into bins and points are colored coded 
by the central oxygen abundance of each bin. Only points with 
measurements of both [S III] and [O III] temperat ures are plotted 
as solid circles. The models of lBinet tg_^t_^yj2Q^). which combine 
the continuous starburst model o f IKewl ev et al.l ()2QQU ^ and the 
shock+precursor model of lAllen et al.l ||2008I) . are plotted as open 
squares. 

The measure of nitrogen across NGC 628 offers fur¬ 
ther insight into the nature of abundance gradients. 
Because nitrogen and oxygen are formed primarily in 
intermediate- and high-mass stars respectively, they 


are returned to the interstellar medium on differe nt 
time scales (e.g. , iGarnettJ Il99fll : Ivan Zee et al.l [T^98allbl : 
iBerg et ^l2012t) . In the upper panel of Figure ITOl we 
plot the CHAOS N/0 gradient, including 11 observa¬ 
tions with T[N II] measurements and three observations 
with T[0 III] measurements from B13. The values from 
the literature are consistent with the trends established 
by the CHAOS observations; the main difference is the 
extension to grea ter radius afforded by the B13 points. 
IBerg et ^ (1201311 used a likelihood ratio F-test to show 
that the N/0 gradient is best described by a bi-modal 
relationship (short dashed line), where decreasing N/0 
flattens out near the isophotal radius (R 25 "^11 kpc,). 

In the lower panel of Figure [10] we illustrate the oxy¬ 
gen abundance relationship across the disk of NGC 628. 
In comparison to the CHAOS data (solid points), di¬ 
rect abundances from B13 and C13 are plotted as open 
squares and open triangles respectively. The reduced 
level of scatter amongst sources when prioritizing T[S III] 
and T[N II] is confirmed when incorporating the ad¬ 
ditional points from the literature. Since the infrared 
fine-structure lines used by C13 are thought to produce 
abundances that are largely independent of electron tem¬ 
perature, the fact that these values lie in unison with 
the CHAOS data offer another conhrmation of the cho¬ 
sen temperature determinations. Further, B13 measured 
T[N II] in 11 H II regions to obtain directly comparable 
abundances, and these values agree with the dispersion in 
the CHAOS sample. In short, this review of the oxygen 
abundance gradient with respect to a variety of temper¬ 
ature determinations supports our decision to conduct a 
homogeneous abundance analysis based on T[S III] and 
T[N II]. 

5.2.2. Comparison with Results from Rosales-Ortega (2011) 
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The most comprehensive emission line abundance 
survey of NGC 628 to date was carried out by 
iRosaies-Orteea et al.l l)201lL hereafter ROll). The ROll 
study performed a nearly complete 2D sampling of the 
disk within the isophotal radius, allowing a variety of 
statistically significant analyses. Complimentary to the 
work presented here, ROll performed an in-depth com¬ 
parison of various strong-line metallicity indicators and 
laid a framework of geometric parameters. 

ROll obtained robust 2D fiber spectroscopic obser¬ 
vations of NGC 628, which were converted to a set of 
96 H II region spectra which are directly comparable to 
the CHAOS observations. With this data set, different 
strong-line abundance estimators were used to explore 
the ef fects of choosin g different calibrations. The ROll 
ff-Te (jPilvuginll2005h analysis of the H II region catalog 
mimics the data collection and analysis presented here, 
and allows the following side-by-side comparison. 

In Figure [TUI we compare N/0 and 0/H radial gra¬ 
dients from the ff-Tg calibration of ROll. The top 
panel of Figure [10] illustrates the general agreement be¬ 
tween the direct abundance and ff-Tg N/0 gradients for 
the optical disk. Within the uncertainty of measure¬ 
ment, the CHAOS best fit (Equation jT] log(N/0) = 
(-0.523 ± 0.034) (-0.077 ± 0.005) x R (dex/kpc)) 

parallels the ROll fit (log(N/0) = (—0.50 ± 0.03) -|- 
(-0.080 ± 0.06) X R dex/kpcJB, in particular, note that 
the gradients are within Icr of one another. Since the 
N/0 ratio has little dependence on temperature, the 
two methods provide nearly identical results except for 
a small average shift. Note the outlined, shaded box dis¬ 
plays the 0.3 dex spread in the measured data as reported 
by ROll. Similarly, the bottom panel of Figure [101 com¬ 
pares these two methods for the 0/H gradient. With 
the addition of data from the literature extending into 
the outer disk (beyond R 25 ), ROll propose a tri-modal 
oxygen abundance gradient: a steeply declining abun¬ 
dance gradient for the majority of the optical disk, with 
flatter features within the inner 2 kpc and in the outer 
disk. The mid 0/H slope (—0.034 ± 0.02 dex/kpc) from 
ROll, extending from ^ 0.1 — 0.7 R/R 25 , is shown as a 
grey shaded box in Figure [TOj whose width represents the 
1(7 error-weighted dispersion in the log(0/H) measure¬ 
ments. For this part of the optical disk, the ROll results 
are in agreement with the slope from the CHAOS obser¬ 
vations (—0.044±0.011 dex/kpc), but are offset lower by 
^0.1 dex. The CHAOS observations do not cover the 
inner or outer regions where ROll find different slopes. 

5.3. Dispersions in Abundances at a given 
Galactocentric Radius 

The dispersion in abundances within the ISM of a 
galaxy is a very important observable, but there are very 
few galaxies with a sufficient number of observations to 
make statistically significant conc lusions. Based on 61 
direct [O III] A4363 measurements. ]Rosolowskv D. Simon! 
(|2008l hereafter RS08) determined a significant disper- 

ROll use a distance of 9.3 Mpc to deproject radial distances, 
in contrast to the 7.2 Mpc assumed here. Since the most reli¬ 
able distance determinations, such as the TRGB distance, have 
not been performed for NGC 628, this choice is simply one of pref¬ 
erence. Here we use the gradient fits in terms of luminosity distance 
from ROll and convert them to match the scale of the CHAOS 
relationships. 


sion in the oxygen abundance gradient of M33 (oi = 0.11 
dex). In this study, a number of parameters, including 
L(Hq;), c(H/ 3), Balmer absorption and emission equiva¬ 
lent widths, and ionization correction factor, were tested 
for correlated against the derived metallicity gradient 
and its residuals, but no trend indicative of systematic er- 
rors was found . The findings of RS08 were challenged by 
iBresolinI (|2011ll using 8 observations of the [O III] A4363 
auroral line a nd the R^s. stro ng-line method. Based on 
this analysis, iBresolin] ^20111 1 measured a much smaller 
dispersion, and speculated that the significant scatter 
seen in the Tg-based abundances of M33 from the RS08 
sample is due to measurement errors in the weak [O III] 
A4363 auroral line. 

Although relative abundances from the CHAOS obser¬ 
vations show very small dispersions at a given galacto¬ 
centric radius in NGC 628, the absolute oxygen abun¬ 
dances show a large and significant dispersion. ROll 
also measured a large dispersion in metallicity, finding 
the scatter to be larger in their fiber-to-fiber sample (rms 
= 0.128 dex) than the H II region catalog (rms = 0.070 
dex) (where several fiber spectra are averaged together). 
This level of spread is consistent with the CHAOS sample 
{at = 0.165 dex, at = 0.159 dex). 

In Figure [TT] we plot the average offset in 0/H abun¬ 
dance from the best fit relationship in radial bins of 0.2 
R/R 25 for both the CHAOS sample and the ROll study. 
For the three bins in common between 0.2 — 0.8 R/R 25 , 
the ROll data has both smaller average departures from 
the measured radial gradient and smaller dispersions 
around these mean offsets. This result further suggests 
that the large dispersion measured in the CHAOS gra¬ 
dient has a physical origin whose effect is removed when 
averaged over entire H II regions. Be that as it may, 
until the source(s) of the discrepant temperature mea¬ 
surements are better understood, at this point, it is best 
to consider the measured dispersion to be an upper limit. 

In addition to temperature discrepancies, we investi¬ 
gate other sources of scatter in oxygen abundance, such 
as the symmetry of the metal abundance distributions, 
but find no obvious relationship with spatial structure 
(differences amongst spiral arms or arm versus inter-arm 
abundances) or ionization fraction. One advantage of the 
ROll study over CHAOS is that it samples a nearly com¬ 
plete coverage of the disk. Since we have shown that the 
two studies find comparable oxygen abundance gradients, 
we can use the ROll disk geometry analysis to comple¬ 
ment our understanding of the CHAOS data. ROll per¬ 
formed two detailed spatial analyse s: ( 1 ) by quadrants 
and (2 ) by spiral arms, as defined bv iKennicutt &: Hodeel 
(|1980ll . ROll found distinct radial trends in the ion¬ 
ization parameter and metallicity for both the geomet¬ 
rical quadrants and the morphological arms, indicating 
that the disk may have evolved somewhat asymmetri¬ 
cally. However, these variations result in 0/H gradients 
which lie within the uncertainties of the overall abun¬ 
dance trend for the disk. 

While the CHAOS and ROll studies suggest the sig¬ 
nificant intrinsic oxygen abundance dispersion found for 
NGC 628 is real, a definitive source for such variation at a 
given radius has not been found. Further, consider i ng the 
the contrasting results of the RS08 and IBresolinI (|201ir ) 
studies for M33, we speculate that these azimuthal varia¬ 
tions are simply not recovered in small samples. For the 




























18 


Berg et al. 

CHAOS Abundance Gradients for NGC 628 



Fig. 10.— N/O and 0/H plotted versus galactocentric radius for NGC 628 and are compared to direct and ff-Te methods in the literature. 
Data from B13 (open squares) and from 013 (open triangles) are plotted, where electron temperatures were determined using [O III] and/or 
[N II] for the former and the infrared fine-structure [O III] 88/i line for the latter. The solid line depicts the least-squares best fits to the 
CHAOS data, as reported in Section [4] whereas the shaded box represents the negative gradient and dispersion measured for the optical 
disk by ROll using the ff-Te method. In the top panel, the CHAOS data intersects the bi-modal N/O relationship observed by B13 (dashed 
line). In the bottom panel, the [N II] based B13 data and the nearly temperature independent measurements of C13 coincide with the 
CHAOS data. 


large samples of direct abundances presented here and 
in RS08, the significant intrinsic dispersions measured 
are due to the electron temperatures. However, the dis¬ 
persion could be due to discrepant temperatures, uncer¬ 
tainty in the measureme nts, or both. Followi ng the line 
of argument laid out in iBinette et al.l l)2012f) and sum¬ 
marized in Section O the primary factor influencing 
temperature discrepancies, and subsequently determined 
abundances, is not yet clear. 

6 . SUMMARY AND CONCLUSIONS 

The CHemical Abundances of Spirals (CHAOS) 
project seeks to establish a broader understanding of the 
chemical evolution of spiral galaxies in general. CHAOS 
harnesses the combined power of the Large Binocu¬ 
lar Telescope (LBT) with the large spectral range and 


sensitivity of the Multi Object Double Spectrographs 
(MODS) to observe “direct” abundances in a large sam¬ 
ple of H II regions in spiral galaxies. In this manner, we 
measure the largest sample of the highest quality spec¬ 
tra to date in NGC 628, with one or more temperature- 
sensitive auroral lines ([o III] A4363, [N II] A5755, and 
[S III] A63I2) being observed at a strength of 3a or greater 
in 45 H II regions. Analysis of this data set provides the 
following insights: 

1) Due to the obvious advantages in observing [O III] 
A4363 over the other auroral lines, it has long been the 
principle nebular temperature diagnostic. Further, the 
scarcity of observations of other lines in the literature 
has postponed any comparative tests of validity. The 
CHAOS sample presented here records multiple auroral 
lines in many of the H II regions observed, allowing a 
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Fig. 11. — Average offset in oxygen abundance from the mea¬ 
sured radial gradient plotted in radial bins of 0.2 dex for each the 
CHAOS sample and the ROll study. Colored solid lines represent 
the average offset in each bin, where the extent of the shaded box 
depicts the dispersion around the mean offset. 

comparison of derived temperatures from multiple emis¬ 
sion lines. We have shown a very tight relationship be- 
tween temperatures based on [S III] A6312 and [N II] 
A5755 indicating that these temperatures are reliable. 
In contrast, temperatures from [O II] AA7320, 7330 and 
[O III] A4363 can often show large discrepancies com¬ 
pared to T[S III] and T[N II]. 

2) These results lead us to adopt a uniform method 
of electron temperature determination in which T[S III] 
and T[N II] are prioritized and used in combination with 
the theoretical relationships of G92. Through this ap¬ 
proach, we were able to minimize the dispersion seen in 
the abundance gradients. 

3) Relative abundances of S/0, Ne/0, and Ar/0, 
which are nearly independent of temperature, are mea¬ 
sured to be constant across the entire optical disk of 
the galaxy. Using weighted least-squares fits, with er¬ 
rors in both the x and y variables, all three show very 
small intrinsic dispersions (0.08, 0.09, and 0.04 respec¬ 
tively). These unvarying values over the range in metal- 
licity are consistent with a universal value for the upper 
mass range of the IMF. 

4) Radial gradients are measured for both N/0 and 
O/H out to the optical radius of NGC 628. The N/0 ra¬ 
tios show a negative gradient, cons istent with the previ¬ 
ous wo rks of iBerg et al.l (|20I3ll and iRosales-Ortega et al.l 
(j20IItl , but do not sample the radial extent necessary to 
assess the presence of a bi-modal relationship. The small 
measured dispersion about the N/0 relationship is at¬ 
tributed to intrinsic dispersion {at = 0.083; at = 0.060). 
We interpret this result as an indication that the ISM in 
NGG 628 is chemically well mixed. 

5) O/H abundances show a moderate negative radial 
gradient of —0.485 ± 0.122 dex/R 25 and a significant dis- 
pe reion of at = 0.165 de x , in li ne with the previous work 
of iRosales-Ortega et al.l (1201111 . We posit that this dis¬ 
persion represents an upper limit to the true dispersion 


in O/H at a given radius. The dispersion in the N/0 
gradient is signihcantly smaller than in the O/H gradi¬ 
ent. Since N/0 is nearly temperature independent, some 
of the scatter seen in the O/H relationship is likely due 
to systematic uncertainties arising from the temperature 
measurements. 

6) Currently, typical analyses of H II region abun¬ 
dances assume a simple symmetric, 3-zone model. The 
electron temperatures determined for these zones show 
large discrepancies, but their physical nature has not 
been resolved. Factors such as shocks, temperature 
inhomogeneities, abundance inhomogeneities, and non- 
Maxwellian electron energy distributions may each play 
a role. Improvements in atomic data, photo- and shock- 
ionization models, and further high-quality, statistically 
significant observations are needed to address this issue. 
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APPENDIX 

OXYGEN ABUNDANCES USING T[0 III] 

We discussed the discrepancies between [S III] and [O III] electron temperatures in Section [5.1.11 and charted them 
in PigureO The uniform analysis approach adopted in this paper prioritizes calculations of the high ionization zone 
temperature from T[S III] over T[0 III]. Here we offer an additional validation of this choice. 

In Eigure[T2 for the 18 regions with [O III] auroral line detections, we plot the oxygen abundance gradient adopting 
[O III] as the standard electron temperature indicator. The solid black line is the error-weighted least-squares gradient 
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Fig. 12.— The oxygen abundance gradient for NGC 628 u sing only [O III] for temperature determinations and the subsequent abundance 
analysis. The best least-squares fit determined in Section 14.3! is compared as a solid line, and the dispersion about this relationship is 
indicated. Regions used in the temperature discrepancy analysis are denoted as colored points, where the highly discrepant points are red 
and non-discrepant points are blue. Log(Ar/0) outliers are outlined by green boxes. 


(Equation [8|) that we determined in Section 14.31 Regions with both [O III] and [S III] temperature determinations 
are color coded according to magnitude of the temperature discrepancy, and log(Ar/0) outliers (as discussed in Ap¬ 
pendix |B]) are outlined by green boxes. Several trends are notable in this plot: (1) the scatter in oxygen abundance 
seems to be systematically shifted to lower abundances, (2) the greatest deviants from the gradient have large tem¬ 
perature discrepancies, and (3) the disp ersion relative to the best fit is large {at = 0.301 dex). Thus, the adoption of 
our analysis approach (see Section l4T|) resulted in a a significant reduction in the scatter {at = 0.165 dex). Because 
Ar and S abundances are sensitive to an intermediate temperature, a propagation of discrepancies ensues if [S III] is 
in disagreement with T[0 III] and T[N II]. Minimizing these temperature inconsistencies removed many of the Ar/0, 
S/0, and 0/H disparities seen in our data. 

INDICATORS OF TEMPERATURE DISCREPANCIES 

[O III] has been and continues to be a useful temperature indicator, but significant deviations have emerged such that 
multiple cases must be considered, and the appropriate modifications applied: (1) In the case that T[0 III]'^T[S III], 
either indicator or an average may be adopted. (2) When T[S III]-T[0 III] discrepancies are large, the multiple-line 
diagnostics presented throughout this paper indicate T[S III] should be used. (3) When T[0 III] is the only measured 
temperature some risk of overestimating the temperature may be introduced, but in a systematic way such that oxygen 
abundance is underestimated. 

In the last case of only one temperature determination, there is no potentially better alternative. However, if multiple 
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Fig. 13.— Temperature sensitive indicators log(Ar“^"*"/S'^'*') and log(77) are plotted for NGC 628. Colored points are used to indicate 
regions where both T[S III] and T[0 III] were measured, where highly discrepant regions are colored red and non-discrepant regions are 
blue. Regions with only one temperature measurement are colored black. Log(Ar/0) outliers are outlined by green boxes. 

auroral line detections are present, a further indication of how to proceed is needed. In this section we present two such 
indicators that appear promising as possible temperature diagnostics. We consider the standard choice of temperature 
indicators: T[0 III] for the high ionization zone, T[S III] for the intermediate ionization zone, and T[N II] for the low 
ionization zone. The following plots are intended as an introduction for analysis to come with the completed CHAOS 
dataset, and so no conclusions are drawn here. 

The ratio of [Ar III] A7136 to [S III] A9069 has been found to be insensitive to ionization s tructure, determined a lmost 
completely by the electron temperature and the ratio of the number of Ar and S ions (|Stevenson et al.l 119931) . and 
therefore offe rs a more controlled lo ok at temperature discrepancies. In Figure [13] we plot Ar++/S+“*' ratio versus the 
77 parameter (|Vilchez fc PagellFl988ll for H II regions in NGC 628 . The rj parameter {rj = log[(0+/0++)/(S+/S++)]) 
represents the hardness of the input ionizing spectrum, and is a sensitive function of stellar effective temperature. This 
parameter pairing allows a more isolated temperature sensitive comparison and the implied result is striking. Most of 
the data follows a a linearly increasing relationship for log(Ar++/S+“'') vs log( 77 ), which is further confirmed by the 
non-discrepant T[S III]-T[0 III] blue points that lie neatly along this trend. Further, most of the highly temperature 
discrepant red points are clear outliers. Some black points are also outliers, but they do not contain both temperature 
measurements required to assess temperature discrepancy. Based on this plot, Ar“''+/S“'"^ is a useful diagnostic to 
identify extremely discrepant values of T[0 III]. 

A second temperature discrepancy indicator of interest can be found in the a-element abundance ratios. In general, 
a-elements are thought to be produced primarily in massive stars and thus returned on similar time scales such that 
their ratio remains constant across oxygen abundance. However, Ar“'"+ and 0^“'' lie in different ionization zones, and 
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Fig. 14.— The Ar/0 ratio as a function of oxygen abundance is plotted for NGC628, where T[0 III] is prioritized for the high ionization 
zone. The expected a-element trend is recovered for the majority of the sample, but the resulting Ar/O value s are higher on average 
(log(Ar/0) = —2.033) than the average value from our standard analysis (log(Ar/0) = —2.219; see Section 14.211 and the scatter greatly 
increased. A number of high Ar/O outliers (log(Ar/0) > —2.0; above dotted line) are present as indicated by green boxes. It is of interest 
that many of these outliers correspond to large discrepancies in T[S Ill]-T[0 III]. 

so if large temperature discrepancies exist, we can expect to see regions with outlying values of log(Ar/0). Figure [14] 
demonstrates the usefulness of the Ar/O ratio as an indicator of temperature discrepancies. Log(Ar/0) is plotted 
versus oxygen abundance, using the standard choice of temperature indicators in the abundance analysis, and for the 
majority of the sample the anticipated constant trend is recovered, although the error-weighted average value is shifted 
up to log(Ar) = —2.03 and scatter is more than doubled {at = 0.186; at = 0.172). Above the dominant horizontal 
distribution of non-discrepant temperature points are a number of points with log(Ar/0) > — 2 , whose separation is 
denoted by a dotted line. Many of the Ar/O outliers, denoted by green boxes, correspond to large measured temperature 
discrepancies, and appear to be systematically offset to higher Ar/O ratios from the main sample. These outliers are 
consistent with the picture that overestimated [O 111 ] temperatures produce underestimated oxygen abundances and 
thus overestimated Ar/O ratios. From these analyses, Ar has emerged as a possibly important gauge of temperature 
validity, and will be revisited in a more thorough analysis in a future paper. 
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